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INTRODUCTION 

Today as petrologists, mineralogosts and crystallographers, we are very 
well equipped in our laboratories. We have refined instruments for chemical 
analysis such as laser probe, electron microprobe and the ion microprobe mass 
analyzer which can produce results within minutes, results which used to take 
days and weeks. Besides x-ray methods, there are several spectroscopic tech- 
iques, such as MOssbauer spectroscopy, which can be used to know the distribution 
of cations over the nonequivalent structural sites in a crystal. The experimental 
methods of synthesizing minerals assemblages in the laboratories have improved 
very much. There are refined techniques to control the fugacities of gases and 
create low to very high P and T conditions in our petrological experimental . 

These technical capabilities are further enchanced by the use of computers 
which can take care of our numerical data and with the help of suitable thermo- 
dynamic theory tell us about the consistency or inconsistency of our experiments. 

As a result of our augmented experimental capability, we have gathered a 
lot of phase equilibria data both from our efforts to synthesize minerals assem- 
blages in the labortory and from our chemical analysis of natural minerals 
assemblages. In order to make a meaningful analysis of both these types of 
data, it is imperative that we have thermodynamic data on the rock-forming phases 
and crystalline solutions. Unfortunately the progress in obtaining these quantilies 
by thermochemical and calorimetric methods has been very slow and it is 
desirable that we consider the possibility of obtaining 3ueh quantilies by other 
methods. 
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We may consider in suitable cases of retrieving thermodynamic data from 
the phase diagrams themselves. These quantilies may be considered significant, 
if they are obtained consistently from different phase diagrams. However in 
the case of crystalline solutions, thermodynamic functions of mixing cannot be 
obtained without the use of certain solution models. Two types of phase diagrams 
may be considered. The first type is the diagram with the solvus or binodal 
curve bounding a binary two-phase region. The analytical methods have been 
discussed by Thompson (1967), Thompson and Waldbaum (1969) and Green (1970). 
These methods, however, are of limited use for rock-forming silicates, since 
most Fe 2+ -Mg 2+ crystalline sc .Unions do not show any solvr.s relationship. The 
second type of phase diagram is the Roozeboom type figure where we plot the 
concentration of a component in one phase against the concentration of the same 
component in the a coexisting phase. Such distribution relationship, based on 
simple ion-exchange reactions, were discussed initially by Ramberg and De Vore 
(1951) followed by Kretz (1959) and Mueller (I960). It is possible to retrieve 
useful thermodynamic information from such distribution data in ion-exchange 

as collected by Nafziger and Muan (1967), Larimer (1968) Medaris (1969) and 

*•> 

Schulien et al, (1970) among others. 

Useful thermodynamic information may also be obtained by considering 
distribution of cations within the crystalline solution. Ghose (1961) found the 
interesting Fe 2+ -Mg distribution in eummingtonite. Since then such Fe 2+ -Mg 
order-disorder has been studied in several silicates by crystallographers. 
Theoretical framework for considering the homogeneous equilibria of intra- 
crystalline cation distribution has been presented in several papers by Mueller 
(1962), Matsui and Banno (1965), Perchuk and Ryabchikov (1968), Thompson 
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(1962), Matsui and Banno (1965), Perchuk and Rvabchikov (1968), Thompson 
(1969) and Griver and Orville (1969),, 

The purpose of this work is to discuss these methods of obtaining thermo- 
dynamic quanlities and discuss some aspects of partitioning of elements in 
coexisting phases. We do this by considering the definition of crystalline solu- 
tions, the definition of components in a silicate minerals and the definition of 
chemicals potentials of these components. We also consider the solution models 
involved. There are examples on calculating thermodynamic functions of mixing 
in the systems CaW0 4 -SrW0 4 , olivine-chloride solution, and orthropyroxene. 

It is hoped that this work will generate enough interest among the experi- 
mentalists to gather useful ion-exchange data on coexisting phases and among 
the crystallogrphers to gather data on site-occupancies in the rock-forming 
silicates. 
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Section 1 


THERMODYNAMIC RELATIONS IN CRYSTALLINE SOLUTIONS 

Thermodynamic relations between the concentration of a component in a 
solution and its chemical potential and other thermodynamic functions of mixing 
are presented here. The details of the simplifying assumptions and the methods 
of statistical thermodynamics have been given by Denbigh (1965), Guggenheim 
(1950, 1967), Prigogine and Defay (1954) among others. Recently Thompson 
(1967) also considered the properties of simple solutions. Besides presenting 
a summary of thermodynamic relations in binary, ternary and multicomponent 
solutions, we shall be particularly concerned with the difficulties encountered 
in their application to silicate minerals. Some of these problems such as the 
choice of a component and definition of its chemical potential in a silicate have 
been discussed by Ramberg (1952a, 1963), Kretz (1961) and Thompson (1969). 

Crystalline Solutions 

Crystalline solutions considered here are the rock forming silicates form- 
ing isomorphous series with one another. Such crystalline solutions have a 
definite structural frame work with generally two or more kinds of nonequivalent 
structural sites. The type of sites and the ions that occupy them vary in different 
crystalline solutions. The overall crystal symmetry of a solution does not change 
as a function of the composition, though certain microscopic details within the 
crystal i. e. the form and size of the individual structural sites may change with 
changing composition. 

We may consider orthopyroxene (Mg, Fe) 2 Si 2 O g as an example. In the 
crystal structure there are single silicate chains parallel to the c -axis held 
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together by the octahedrally coordinated Mg 2+ and Fe 2+ . There are two kinds 
of structurally nonequivalent sites Ml and M2 occupied by Mg 2 * and Fe 2+ . The 
Ml octahedral space is nearly regular polyhedral but the M2 space is quite 
distorted. As a result of varying Mg and Fe 2+ in the composition of the crystal, 
the general symmetry of the crystal does not change but there are distinct 
changes in Ml and M2 polyhedra. The former becomes more regular and the 
latter more distorted with increasing Fe/Mg ratio. Such microscopic changes 
at the structural sites within the same crystal framework may be regarded as 
continuous and the resulting energy changes a consequence of the mixing or 
solution of the species to form a crystalline solution. 

Choice of a Chemical Component 

The definition of a component in a mineral is not unique. The components 
in orthopyroxene may be considered to be the molecules MgSiO and Fe Si O , 
or MgO, FeO, and Si0 2 , or the cations Mg 2+ , Fe 2+ , Si 4+ and O 2 ". In petrologi- 
cal studies, the choice of a component is determined by known or postulated 
chemical reactions involving a mineral. In such studies, the choice of compo- 
nents such as FeSi0 3 or FeO is convenient even though there are no discrete 
units of this kind in the orthopyroxene crystal structure. However, while con- 
sidering the thermodynamic properties of 'silicate crystalline solutions, it is 
only realistic to consider the ions as the components (see Bradley, 1962). In- 
deed it can be noted that if the substitution of the cation Fe 2+ by Mg 2+ in orthro- 
pyroxene does not induce any changes in the silicate framework or if there are 
any slight changes, they are directly a function of the changing Fe/Mg ratio, the 
alternative methods of defining FeSiO a or Fe 2+ as a component are equivalent 
(see also Saxena and Ghose, 1971). 
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Chemical Potential and Activity of a Component in a Mineral 
A solution is ideal if the chemical potential of every component is a linear 


function of the logarithm of its mole fraction according to the relation 


fi ^ = p 1 + R T 1 n x t 


( 1 . 1 ) 


where /a. is the chemical potential of i in a solution and p. 1 is the chemical 
potential of pure i. p 1 is a function of P and T only. In a binary solution a , 
whose composition is (A, B) M where M may represent the anion group or the 
silicate framework, and A and B the cations which substitute for each other, 
there is a choice between adopting the cations A and B as components or the 
end member molecules AM and BM. As noted before, under certain conditions, 
the mole fractions may be calculated as 


x ^ - 


A A + B 


or x * = 


AM 


** AM + BM 


and these could be considered equivalent. We may write for chemical potentials 


p? = P AM + R f 1 nx * 


( 1 . 2 ) 


or 


^am = /**“ 


AM 


( 1 . 3 ) 


where p*** and are chemical potentials of A and AM in a standard state. The 

A 

standard state AM is well defined but the standard state with reference to cation 
A needs definition. In orthopyroxene, this is tike referring to the chemical 
potential of Mg 2+ in pure (Mg, Mg) Si 2 O e . The Gibbs free energy for the pure 
end member MgSi0 3 is defined and measurable experimentally but the meaning 
of free energy of Mg 2+ in pure enstatite is little understood and experimental 
methods remain to be developed for its measurement. 
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However in theoretical discussion, where we are not concerned with the 
measured values of the potentials , the definition of chemical potential of a cation 
in a crystalline solution is not only permissible but also useful. Kretz (1961) 
defines the chemical potential of Mg in orthopyroxene as 



where n is the number of cations in the formula. 


(1.4) 


In many crystalline solutions, when their compositions are expressed in the 
simplest form, there are two or more cations in one mole. Examples are oli- 
vine (Fe, Mg) 2 Sib 4 and garnet (Fe, Mg) 3 Al 2 Si 3 0 12 . The chemical potential 
of a component using 'molecular' model is expressed as: 


Fe„SiO. 

^Fe 2 sio 4 + R T In * Fe2 s i04 


01 


(1.5) 


Fe,Al 9 Si,0, „ 

/Li H* -= n 3 12 + ft T In *_ A1 

Fe 3 Al 2 Sl 3°12 Fe 3 Al 2 Sl 3°12 


( 1 . 6 ) 


where the first /x is the chemical potential of the end member in the solution, 
the second /x is the chemical potential of the pure end member. If we choose to 
consider the cation Fe 2+ as component, we have 


=/4e 2Sl ° 4 + 2RT inxi 1 


(1.7) 


Fe 


and 


Fe,Al,,Si ,0, ,, 

=^Fe 3 12 + 3ft T In *i? r 


( 1 . 8 ) 


Fe 


where the first fi is the chemical potential of Fe 2+ in the crystalline solution 
and second ^ is the chemical potential of Fe 2+ in the pure end member. The 
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mole .fractions * are the same quantities in both the ’molecular’ or ’ionic’ 
models. It may be desirable to consider the chemical formula on one cation 
basis i. e. we consider olivine as (Fe, Mg) Si O and garnet as (Mg Fe) AU,,. 
SiO.. In these cases we may write 


_ol „ _1 Fe 2 Sl0 4 

'Fe 2 ’ Fe 


f R T In x. 


(1.9a) 


_ 1 FCjSiO,! 


FeSi o,s°? 2 Fe 


+ R T lnx°{ 


(1.9b) 


gar _ 1 Fe 3 A1 2 Si 3°12 

'Fe “ 3 ^Fe 


+ r t in 


(1.10a) 


trnr 1 Fe , Al ,, Si ,0 , ,, 

^FeAl 2 / 3 Si0 4 = 3^Fe 3 2 3l *+ R T In 


(1.10b) 


The usefulness of the above relations is mentioned later in connection with 
the composition of coexisting minerals. 


We shall be mostly concerned with the activities of the components. For 
a binary ideal solution the activity is equal to its mole fraction. In olivine the 
activity of fayalite (fa) molecule is 


“fa “ 


or for Fe 2 , 


“fe = (*Fe) 2 


Similarly for garnet we may have 


a . rr x *7‘ 
aim aim 

■>„. = (*m 3 
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It is desirable fco consider many reactions, particularly the ion-exchange 
reaction, on a one cation basis i. e. consider olivine as (Fe Mg) Si 0 s 0 2 etc. 
Activity of a cation is then equal to its mole fraction. It is necessary to specify 
that although in this situation x (mole fraction Fe 2+ /Fe 2+ + Mg) is numerically 
the same as * fB (percent of fayalite), the activities are different. * Fe is equal 
to (a pe ) 1/2 in the ’ionic model' but * fn is equal to a fa i. e. activity of fayalite 
in the solution. 


Non-Ideal Binary Solutions 

The relation between the chemical potential of a component i and its activity 
in a solution is given by 


(l i +nlna. 


( 1 . 15 ) 


ideal solution is the limiting case when a. is equal to the mole fraction In 

all other cases, one may express the relation between a L and* i as: 

„ ( 1 . 16 ) 
-i =fi x i 

where f ■ 1 is the activity coefficient of the component i in the solution. 


The free energy of mixing 6' M for a binary solution (A, B) M is given by 

G U= X A RTina A^B RT ^ a B (1 * 1?) 

= RT(x a In x A + x B in * B ) + R T {x A In f A + x B In f B ) 

= C IM + G EM 

The first term G is the ideal free energy of mixing and the second term 
C „ is the excess free energy of mixing arising due to the non-ideality of the 

EM 

system. G m is one of the functions of mixing termed 'excess functions'. Details 
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of the excess functions of mixing may be found in Prigogine and Defay (1954) 
and Thompson (1967). 


Regular Solution Model 

We are often required to calculate the thermodynamic properties of the 
crystalline solutions from compositions of mineral assemblages either obtained 
in experiments or occurring in rocks. In many situations the use of certain 
models for the activity composition relationship helps to assess such properties 
closely. Guggenheim’s (1952) regular solution model is next in simplicity to 
ideal solution model. 


Zeroth approximation 

The excess free energy of mixing G m in a regular solution with the zeroth 
approximation i. e. approximation of complete disorder, is given by 


G = x x W* 

EM A B 


(1.18) 


where A and B are components of a solution (A, B) M and w' is equal to Nw . N 
is Avogadro’s number. W' is often referred to as the ’interchange energy'. We 
shall be very much concerned with regular solutions in this work. Therefore 
we shall briefly discuss the parameter w. A simplified account of this parameter 
is presented by Denbigh (1965). It is assumed that the cations A and B are of 
roughly same size and can be interchanged between lattice sites without change 
of lattice structure and without change in the lattice vibrations. There is a 
certain interaction between A and B given by the energy w which is given by 


W = C 2 "ab “ W aa " tt 'bb ) 


(1.19) 
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where w aa is the increase in potential energy when a pair of A ions are brought 
together from infinite distance to their equilibrium separation in the solution. 
Similarly w ab and u^ b are pair potentials for A-B and B-B pairs. In sp: _ of 
the interaction energy, it is assumed that the mixing of A and B is random. 
This means that the entropy of mixing is the same as that for an ideal solution 
and deviations are expressed entirely in terms of the heat of mixing. 


The thermodynamic equations for the regular solution model with zeroth 
approximation are: 


^EM = ^EM “ *A X B y 


S EM = 0 


( 1 . 20 ) 

( 1 . 21 ) 


The ’interchange energy’ W is independent of p and T. As the excess entropy 
of mixing is zero, according to this model, the predictions on the excess free 
energy of mixing and the heat of mixing (which may often be different from the 
G m ) are not satisfactory. 


Simple mixture model 

In the regular solution model W' is supposed to be independent of temperature 
and pressure. In Guggenheim's (1967) latest version of the lattice theory W 
may be treated as an adjustable constant required to fit the experimental data 
to the model. Such an energy parameter with a symbol W may be called 'a 
cooperative free energy’. 2 W is in a sense the free energy increase in the 
whole system when an AA pair and a BB pair are converted into two AB pairs. 

It is expected that if W is fitted to the free energy data at each temperature, the 
large errors usually found in the predictions of G m and Ii m with composition 
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may be afc least partly eliminated. For a random mixing approximation, the 
various excess functions are given by 


C EM = *A *B W 

e - 3 IF 

"* EM “ X A X B g-yf 

,, /. r „ b ir 

,y EM = *A *B ( 1 ' T 


The activity coefficient is related to the mole fraction as: 


In/. = — (.%- ) 2 


( 1 . 22 ) 

(1.23) 

(1.24) 


(1.25) 


Quasi-chemical model 

The main assumptions required for this model are similar to those of the 
regular solution model in the preceding sections namely a) only configurational 
partition function of the solution is contributing to the thermodynamic excess 
functions and b) the intermolecular forces are central and short ranged and 
therefore the internal energy at 0° K may be obtained by an addition of the pair- 
potentials. The assumption of complete randomness is not required here. There- 
fore any differences found in the calculated values of the excess functions of 
mixing by the zeroth approximation and by the quasi-chemical approximation is 
due to the effect of ordering considered in the latter. 


In binary solutions where the two components A and B are of similar size, 


the activity coefficients are given by the equations: 



73+.1-2*b\ 2 

*a (.3 + 1) J 


(1.26) 
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B 


R - 1 + 2 *, 

* A c/s + *r 


(1.27) 


where z is the coordination number and P is given by 


^{1 + 4 , a ( e 2W/ zRT - l )} 1 ' 2 


(1.28) 


/3 is unity for a perfectly random mixture. /3 > 1 indicates a tendency for clus- 
tering and P < 1 indicates a trend for compound formation. 

The excess free energy of mixing and the heat of mixing are given by 


r 1 P T J , + 1 - 2 * B ft - 1 + 2 * B 

G-... =_ z /? / <.r. In + In 

EM 2 ] A x. (P +1) B x (P + 1 ) 


(1.29) 


H em - ^ + 


if - T —"l 
3 + 1 Ab \ 3 7 J 


(1.30) 


The various equations of the quasi-chemical approximation may be expanded 
as power series in %W /zRT: 


6 em ' w 


X -V il - — 
A B \ 2 


1 / 2 W \ 1/2 WY , . o 

2 rn) ** *» ‘ e tarr ^ * b <*« - *b > 2 + 


} 


(1.31) 


/ A - Wx*il +i (1 -30 + i(AE_| x. 


2 Vz R Ti 


6 \z R T j 


(1 - 11 x B + 28 - 20 x*) + . . . 


B' 


} 


Similarly / B may be obtained by replacing A by B in (1.32). 


(1.32) 
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For molecules which are not very similar in size, we find the necessity to 
include a ’contact factor' (see Guggenheim, 1952 p. 186) in the above equations to 
take into account the size differences. The contact factors may be found 
roughly proportional to the molar volumes or ionic radii. The activity coeffi- 
cients are given by: 


/a 


1 + 


</>B O 3 " 1) 

(£ + 1) 


ZC V 2 


(1.33) 



</> A C/3 - 1) 
0 B C/3 + 1) 


zq 


b/2 


(1.34) 


where q A and g B are contact factors related to the contact fractions <P A and 0 B 
and the mole fractions x A and x by the following equations 


*A 


*A ?A 


*A H A + *B ?B 


c p ~ 


X B ? B 


X A 9 a + x h9b 


(1.35) 


For more details on the derivation and significance of the constants q and g 
and the fractions <f>. and <£_ , reference may be made to Guggenheim (1952, p. 186) 

A D 

and King (1969), p. 488). (5 in (1.33) and (.1.34) is obtained by replacing x. and 

A 

x B by <p A and <£ b respectively in (1.28). 

The other excess functions are given by 
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(1.37) 


«„ -- llihl&L — fr-ril 

<**»* + S’b* <■" + U V : ' 7 


.S’ can be obtained by the standard equation 

EM 


f, EM ^EM ” 1 ^EM - 


Some Other Methods for Representing the Activity-Composition Relations 
In the case of a binary system of components A and B, the Gibbs -Duhem 
equations is 

*A din / A + * B d ln/ B = 0 (1 ’ 38) 

The changes d Inf A and d lnf Q when due to composition change d* A at constant 
temperature we may write 


3 In f 


Bln f 


x , 


B x 


+ x 


B 


B 


= 0 


(1.39) 


A /T 


x 


a y t 


A solution to the above equation was proposed by Margules in the form of a 


power series 


ln - “» + ‘a *8 + °A + rf A *B + • ' ’ 

ln = °B *A + “b X l + C B *1 + d B X i + • • • 


(1.40) 

(1.41) 


when the series is terminated at * 3 , the following relations exist between the 
coefficients 

£3 

• A = = ° (1 ' 42) 
b B ~ b A + 3 G A //2 + 2 d A + . 

C B ~ " (c A + 8 d A /3 + • • 0 
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Using the above relations, Carlson and Colburn (1942) expressed the activity 
coefficients by the equations 


log f A = (2 H - ,1) (1 - * A ) 2 + 2 (A - B) (1 - x A y 
log f B = (2 A - B) v 2 (B - A) x 3 


(1.43) 


(1.44) 


Relations similar to these have been used by Thompson (1967) and Thompson 
and Waldbaum (1968, 1969). 


Another two-constant equation is due to van Laar. The equation resulted 
from a theory based on the van der Waals equation of state. This theory is 
probably incorrect but van Laar's equation continues to be useful for represent- 
ing the activity-composition relation. This equation is 



similarly for the other component 


(1.45) 


io g/ B 


B 



(1.46) 


For many chemical systems the van Laar’s equation provides a better repre- 
sentation of the data than is given by the Mar gules two-constant equation. The 
relative merits of these two equations were discussed by Carlson and Colburn. 

It may be finally remarked that a power series ejqmnsion as (1.47 see later) for 
the excess free energy is now widely preferred and we shall, therefore, use 
only such expressions and not the equations mentioned in this section. Expressing 
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the excess free energy as a power series is a means of giving empirical de- 
scription to deviations from ideality which is a better alternative to the power 
series expansions referring to individual activity coefficients mentioned above. 
0 m expressed as a power series can be related more conveniently to the other 
global properties of the mixture, such as the heat and volume change of mixing 
than can the individual activity coefficients which represent the deviations divided 
up, as it were, among the components. 

f 

General Relations for Binary, Ternary and Quaternary 
Non-ideal Crystalline Solutions 

Excess functions in non-ideal solutions may conveniently be expressed by 
a power series in the mole fraction. Guggenheim (1937) suggested that C EM 
can be expressed as a polynomial in.%- as 

C EM = X A X B H) + A 1 O'A " *B> f *2 (*A ~ + * * *■ > 

where A Q , ,4 and/l 2 etc. are constants. When odd terms in (1.47) vanish, the 
solution becomes symmetric. If and other higher terms are also zero, we 
have the simple mixture model with /t as the energy constant IP in (1.22). The 
expressions for the activity coefficients are obtained from 


R T in f . = + x 


3 a 


EM 


A EM B B 


A (1.48) 

= X l [A Q + A 1 C 3 X A “ * B ) + A 2 ( X A ~ X B ) (5 x a ~ * B ) + 


RTl n /= - x 


3 a 


EM 


B — EM -A 


B 


X A A 1 ( 3 X B ~ X a) + A 2 ( x B - x A ) (5 x B - x^) + 


(1.49) 
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Equations for other excess functions of mixing may be derived from (1.47) 


as: 


- S 


EM 


X A *B 


U 

B — — + 
B T 


fdA,\ 

pT„ 


( B A, 

( x a - * B ) ) (*a - + 


(1.50) 


!L.. — x . X- 
EM A B 


+ ^ 


A o -T 


2 - 1 . 


P A c\ f f d 
\”7 + F ” T \3t7j <a: * ' *»> 


~ 1 V BT 


'<*a - *b) 2 + 


(1.51) 


The expression for the excess free energy of mixing in a ternary system 
according to Redlich-Kister equation (King, 1969) is 

G EM ~ X A X B lt AB * C AB ( X A " *b) + ^AB ( X A ~ X B^ + * '* 

+ X A X C ^AC + C AC ( X A “ + ^AC ( X A ~ X c ) 2 + ' ‘ 


+ *B X C ^BC + C BC Ofi ** + ^BC ( A ‘ B “ * C ) 2 + * ’ 

+ * A *B *C C F + H A X A+ H B * B ) 


(1.52) 


where&.c and d represent A Q , and4 2 of equation (1.47). F ,// A and// B are 
ternary constants and must be evaluated from ternary data. Eor silicate crys- 
talline solutions such data are rare. Therefore we may consider that only the 
pairwise interactions are important and adopt the following empirical expression 
for for a system containing any number of components: 


g em “ 


= / , x i *j f b ij +c ij ( x i -*i) + d ij ( x i * x i) 2 + 


( 1 . 53 ) 


AJ 


where i 4 j. 
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Section 2 


THERMODYNAMIC STABILITY OF A SOLUTION 
Instrinsic and Extrinsic Stability 

A crystalline solution when ideal adds a certain amount of negative free 
energy of mixing to free energy of the system. With increasing positive devia- 
tions from ideality this contribution becomes less and less. Below a certain 
critical temperature of unmixing, the solution unmixes to form two or more 
solutions. These energetic changes obviously affect the stability of the entire 
system of mineral assemblage. This instability of a crystalline solution, which 
arises due to the positive excess free energy of mixing, may be termed intrinsic 
instability (see, Mueller, 1964). Ideal solutions are always intrinsically stable. 
A crystalline solution may also become instable if the physical-chemical con- 
ditions change in such a way that certain reaction products become a lower 
free energy assemblage than the crystalline solution. This instability could be 
considered as extrinsic. A solution may be both intrinsically and extrinsically 
instable. The division is essentially artificial. It, however, helps to understand 
and describe certain petrological reactions as shown by Mueller (1964). 

Olivine (Fe, Mg) 2 Si0 4 and pyroxene (Fe, Mg) Si0 3 may be considered as 
ideal binary solutions at high temperature 1100 °C). In spite of their ideal 

character, orthopyroxenes with more than 55 mole percent of ferrosilite were 
found unstable at liquidus temperatures by Bowen and Schairer (1935). The iron 
rich pyroxene is unstable due to the instability of ferrosilite relative to fayalite 
and quartz. This is the extrinsic instability. 


At low temperatures (~ 600°C), the situation is little different. Orthopyroxene 
is somewhat non-ideal and high values of 6 EM are associated with high ferrosilite 
content of the solution. The extrinsic instability of the solution relative to oli- 
vine and quartz is less since iron-rich pyroxenes (~ 86 percent FeSi0 3 ) are 
stable in metamorphic rocks. The instability of pyroxenes with higher ferro- 
silite in metamorphic rocks may be both due to the extrinsic and intrinsic in- 
stability of the orthopyroxene solution. 


'Critical Mixing 


General conditions: The conditions for critical mixing in terms of free energy 
of mixing and the mole fraction * are 


B 2 C m /B = 0 


(2.1) 


B 3 G m /B * 3 = 0 


( 2 . 2 ) 


These may be expressed in terms of excess free energy of mixing C EM as 

B 2 C em /B * 2 = - R T/x (1 - x) 

B 3 C £m /B * 3 T (2 * - !)/.%• 2 (1 - x) 2 


(2.3) 

(2.4) 


Si mple mixture: For a simple mixture we have 

g em = * a -*) w r=r( t, p) 

By successive differentiation of (2,5), we have 

B 2 Cem /b ^2 = „ 2 r 
B 3 C em /B * 3 = 0 

By substituting (2.3) and (2.4) into (2.6) and (2.7) respectively, we have 


(2.5) 


( 2 . 6 ) 

(2.7) 
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- 2 ir = - R T/x (X - x) 


( 2 . 8 ) 


0 = R 7 (2 X - 1) /x~ (1 - x) 2 


(2.9) 


These give us the critical composition as x = 0.5 and 2 R 7 - If . 

General non-ideal solution: For a binary solution which is not a symmetric 
solution, the excess free energy of mixing is given by 

Cq, = *(»- .V) u 0 , A, (1 - 2 x) A 2 (1 - 2 V) 2 ... (1 - 47) 

Successive differentiation of (1.47) with respect to a gives 


ri 2 a 


EM 


a x 2 


- 2 ,i Q - 6 A t (1-2 V) - 1 2 iio - 48 V (1 - a).- 


( 2 . 10 ) 
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; 48 l 2 (T - 2 x) 


( 2 . 11 ) 


Substitution of (2.3) and (2.4) into (2.10) and (2.11) gives equations which are 
transcendental and cannot be solved without a computer program using an iter- 
ation method. 


Formation of miscibility gaps in a ternary simple mixture: Let us consider a 
ternary simple mixture with components 1, 2 and 3. The W for the three binary 
systems are IF o , If and lf / . The chemical potentials of the components in 
the solution are given by: 

- pO ( T , P) + R T In x 1 + R T In f v etc. , ' 2 - 12 - 

where R Tlnf may be expanded in terms of a- and IF as follows: 
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(2.13) 


R 1 In f 1 ~ (x 2 ) 2 lf 12 4- (* 3 ) 2 * l3 + x j x 3 (V 12 - |f' 23 + H' 13 ) 

RTinf 2 = (x 3 ) 2 W 23 + (xj* r 12 + *3 (r 23 - ip 13 + r 12 ) 

/? 7 In f 3 = (^j) 2 ff 13 + (* 2 ) 2 ^23 + x 2 (^13 “ ^12 + ^2 3) 

At equilibrium in the two separated coexisting phases a and ft we have 


Mi (. x 2 1 * 3 . T) -f4 (x\ t x\, T) = 0 


(2.14) 


and similarly for u 2 and u 3 . Substituting f'2.12) and (2.13) in (2.14) and rearrang- 
ing we have (see Kaufman & Bernstein, 1970, p. 226) 

R T In + IP 12 (xf - xf) + W l3 (xf - xf) 

r A IT (*§ x\ - x* *J) = 0 


« r ^ c*|/*S) + jp 12 [(1 - 4) 2 - (1 - xj) 2 ] + IP 13 (*f - *J ) 

- A IF (1 - *|) - (1 - *J)] = 0 

/? 7’ in + ir 12 (*f - (*j) 2 ) + r 13 [(1 - *f) 2 - a - <) 2 ] 

-A W [x\ (1.- X P) - (1 - *“)] = 0 

where A IF is IF = IF 12 + - JF 23 


(2.15) 


With the help of equations (2,15) compositions of coexisting phases may be 
calculated and the miscibility gap may be plotted on a ternary diagram. How- 
ever first, we require the compositions of the coexisting phases on three binary 
edges. 

In a binary solution, the miscibility gap can be calculated by finding the 
composition of the coexisting phases which together represent the minimum 
free energy of the system. This may be done by the "Tangent method" graphic- 
ally i. e. by drawing a tangent through the two points representing the two 
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- m 


minima in the plot oi' the free energy of mixing against composition alternatively 
we may consider the relations 


• U 1 ~^ l J u 2~^2 

For the binary regular solution, there is a symmetric miscibility gap and there- 
fore 


= 1, x@ + - 1 and also = x@ 


and - xP . We may therefore write 


R T In (1 - xj + x\W ~RT In (1 - * 2 ) + x\ W 
R T in x x + (1 - x t ) 2 W = R T In * 2 + (1 - * 2 ) 2 I V 


Substituting 


'2 = 1 “*1 


in (2.16) we have 


(2.16) 

(2.17) 


If' 


n (i-2 


ln< 


c 1 -*i)l 


x . 


(2.18) 


Equation (2.18) may be solved by an iteration method to find the miscibility gaps 
on the binary edges in a triangular diagram. 


A computer program to solve the equations (2.15) numerically and 
method to form a. miscibility gap have been presented by Kaufman and Bernstein 
(1970). Some examples to illustrate the possible solutions of certain mineral- 
ogical problems are presented elsewhere. 
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Section 3 

COMPOSITION OP COEXISTING PHASES 


Ideal Solution Model 

Distribution of a component between two ideal binary crystalline solutions 
Although there are no strictly binary silicates, certain minerals such as 
orthopyroxene and olivine may be assumed to be quasi-binary. In view of the 
similarity of Pe 2+ and Mg 2+ in ionic charge and size, we may assume that 
olivine and orthopyroxene are binary ideal solutions. This assumption is re- 
examined later. 


Ramberg and DeVore (1951) considered the following ion-exchange equilib- 
rium between olivine and pyroxene: 

Mg Si 0 3 +4 Fe 2 Si 0 4 v^Fe Si C> 3 + *Mg 2 Si 0 4 (3a > 

en fa fs f o 

The equilibrium constant for the above reaction at a certain P and/ is 


K 


(i - 4‘j 


Fe ' 


3a 


(1 - X<£*) 


Fe 


(3.1) 


The equilibrium constant K is a function of P and T only. In the present case 
K. , however, is not found to be constant except at high temperatures (see Olsen 

o fl 

and Bunch, 1970). 


It may be noted that (3a) is written on one cation exchange basis. We may 
also write 

2 Mg Si 0 3 + Fe it Si 0^ 2 Fe Si 0 3 + Mg 2 Si 0 4 ( 3b ) 
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The equilibrium, constant for the above reaction is 


3b " (i - *?r ) 2 

A Roozeboom figure with such K values has been presented by Kern and Weisbrod 
(1967), p. 224). It is known empirically from the distribution data in several 
mineral assemblages, that equilibrium constants or distribution coefficients 
such as K b are very cumbersome to handle and inconsistent with petrological 
observations. One may, therefore, prefer to use the distribution data on one 
cation exchange basis. It is obvious that in actual calculations of the energy 
values , we shall be required to adjust for the activity-composition relations 
(such as 1.7, 1.8, 1.12 and 1.14) as discussed before. 

Generally olivine and pyroxene coexist with several other minerals of 
fixed or variable composition. If there is no significant change in the concen- 
tration of the minor components changing the binary character of the two 
minerals, & 3a is not a function of any changes in the number or proportion or 
composition of other coexisting phases. This is generally true about equilibrium 
constants in other systems also. At a certain P and / the stability of the as- 
semblage olivine-pyroxene is a function of the presence or absence of quartz 
but the value of /C, itself is not affected. 

Kretz (1959) used Roozeboom plots extensively to show the orderly distribu- 
tion of cation between coexisting silicate minerals in rocks. If chemical equil- 
ibrium is closely approached in the distribution of a component between two 
binary solutions at a certain P and the distribution isotherm is a smooth 
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curve and if at the same time both the solutions are ideal, it will be of the 
form shown in Figure 3.1. 



v 

& 






Coexisting ternary ideal solutions 

Consider two coexisting ternary phases a, and /3 with the formulae (A, B, C) 
M and (A, B, C) N. The pure components are AM, BM and CM in a and similarly 
AN, BN and CN in/3. The chemical potentials of the components in a and ft 
are 

(3>3) 

+ in*& <3 ' 4) 

ans similarly others. The potentials of the pure components fx m etc. are 
functions of Pand T only. We may also substitute x* for and for x$ 

J J A AM A AN 

without altering the results (see discussion before). 


The distribution of A between a and P> may be represented by the ion ex- 
changes 


Aa + B /3^^B a + A /S 
A a + C /3^^ C a + A /3 

We may write the equilibrium constants as 


(3c) 

(3d) 


K 


3C 


X K 

B A 
A B 


and = C A 


3d 


*A X C 


where X A is A/A + B + C etc. 

Both K 3C and K 3d will be constants for all ratios of A: B: C. A plot of * k against 
x^ k will produce a symmetric ideal distribution curve. 
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Non-ideal Solutions 


Distribution of a compo nent between two "simple mixtur es" 

For the ion -exchange equation 

Aa + B/6 = A/3 + Ba 

we have at equilibrium 

=l4 

If (A, R) a and (A, B ) .6 are ’simple mixtures', we have 

^ =M Ba + R T In (1 - *J) + W (*J) 2 

and similarly others. 

Substituting such values (3.6) in (3.5) and rearranging we obtain 


(3e) 


(3.5) 


(3.6) 


In 


(! -*£) 



(1 




- AC« 
/? 7 1 


( 3 . 7 ) 


where 


A C J = jj? a + i±W _ - /'A a . 


Or we may write 


. IF a 

lnX3. = ln^-^_ 


(1-2 *J) 


r/ 3 


(1 - 2 


(3.8) 


where & 3e is exp(-AG°//?F)and is the distribution coefficient. 

If we can obtain a good least squares fit for the distribution data by using 
(3.8), it may also be found that both the minerals are close to being 'simple 
mixtures’. 
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With increasing non-ideality in one or both of the minerals, the distribution 
isotherms may attain different forms (see Mueller (1964). Figure 3.2 shows an 
example where one of the minerals a is ideal and 0 non- ideal. \yP is assumed 
to vary linearly with 1/7’. The values of if'/ 3 and K at 673 and 1173°K are 
2.75, 1.603, 0.77 and 1.518 respectively (see Saxena, 1969a). The forms of the 
distribution isotherms are very different from the symmetric ideal curves. 


Coexisting regular ternary solutions 

The composition of two coexisting phases which obey the same equation of 
state is considered here as an example. These phases are products of unmixing 
in a ternary solution (A, S, M. We assume that IP , IP and IP are 1500, 
7000 and 9000 cal/mole respectively and IPs are not a function of p,f and com- 
position (regular solution). Figure (3.3) shows the miscibility gap in the system 
and the tie lines for the coexisting phases. Let the phase rich in C be denoted 
by a and the phase poor in C by 0. For the chemical potentials we have 

/Ij = u m + R T In / 2 x\ (3,3) 

=M am + R T in fPxg etc., (3,4 * 

We may consider any one of the following ion-exchanges between a and 0 


A 0 -!- B a^A cl + B 0 
Ba + C/3 = B/3 + Ca 
Aa + C^CatA/S 


(3f) 


The equilibrium, constant for (3f) is 


'f a fp 

X A X B > A 1 B 
• . _ 

*A f'it l 


= K 


(3.9) 


3 f 
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In this particular case since a and /? obey the same equation of state we have 
AC 0 = 0 and IC { is unity. In other cases where a and /3 are minerals with dif- 
ferent crystal structures, the equilibrium constant is not unity. The f terms 
in equation (3.9) are functions of P, T and the A: B: C ratio and therefore K d 
(X^ xf /X“ Xg) also changes with P , T and A:B, B:C and A:C ratios. 

Let A:B:C change systematically as listed in Table 3.1. We find that a plot 
of x* against x@ where * is either A/A +B or A/A +B +C shows a smooth distri- 
bution curve (Figure 3.4). The form of the curve, however, is markedly different 
from the ideal distribution curve. 

The activity coefficients are given by 

R T ln U = (*B> 2 + (*■/ W ‘>C + *c ^ + “'ic) 

R T 1" f B = (*c) 2 "bC + C-A> 2 n* + *c *A ("Sc - «Sc + "Sb> (2,13) 

R T ln fc = (*a ) 2 *ac + (^b ) 2 "Sc + *A *B <"ac - »'ab + iTSc) 

where x A = A /A + B + C etc. It may be checked that substitution of / values 
into (3.9) gives the equilibrium constant as unity. 

Distribution of a Cation Between Two or More Multicomponent Minerals 

Many rock-forming minerals are complex multicomponent crystalline 
solutions. The distribution of cations in two or more coexisting minerals in 
natural assemblages may still yield certain valuable informations. The method 
to be followed in such cases has been discussed by Kretz (1959). In silicates 
there are at least two types of coordination for the cations. Si 4+ , Al 3+ , Fe 3+ 
and less commonly Ti 4+ are in tetrahedral coordination. Fe 2+ , Mg 2+ , Fe 3+ , 

Al 3+ , Mn 2 *' and Ti 4+ are found in the octahedral coordination. Such differently 
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coordinated ions may be regarded as forming sub- mixtures. The distribution 
of Fe 2+ or Mg 2+ or any other octahedraily coordinated ion may be examined in 
two or more such sub-mixtures forming parts of different minerals. It should 
be noted, however, that the chemical potentials of a cation in octahedral coordi- 
nation may also be a function of any chemical variation in the concentrations 
of the tetrahedrally coordinated ions. Such information can be usually obtained 
before hand by considering the chemical composition of individual minerals. 

For example the positive correlation between the concentration of tetrahedrally 
coordinated Al 3+ in amphiboles and biotite with the Fe 2+ /Mg ratio in the mineral 
is now well known (Ramberg, 1952b, Saxena, 1968a). 

It may be argued that the study of the distribution of a component between 
only two of the coexisting minerals which are quasi-binary solutions out of an 
entire assemblage of five or six minerals cannot be useful. In other words, the 
presence or absence of a third or fourth mineral in the assemblage would affect 
the distribution coefficient. This is not generally true. The distribution co- 
efficient changes only when the presence or absence of a third mineral is as- 
sociated with a significant change in the concentration of one or more elements 
in one or both of the coexisting minerals. TiO s is only sparingly soluble in 
olivine and orthopyroxene. In such a case the chemical potential of Ti0 2 may 
increase or decrease in the rock and rutile may be added or removed from the 
assemblage, K. for the distribution of Fe 2+ and Mg does not change. However 
if the change in^. changes the concentration of TiO significantly in one of 

i 10 2 2 

the two coexisting minerals, may also change. Thus it is only meaningful to 
consider the concentrations of all the components in the two minerals and not 
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the presence or absence of another phase or the change in the bulk composition 
of the rock. 

One of the important results of the study of cation partitioning is the recog- 
nition of how closely chemical equilibrium may be approached in the rocks. 
Irrespective of whether the minerals are ideal or not, the distribution of a com- 
ponent between two coexisting binary phases at a certain P and T will be repre- 
sented by a smooth distribution curve provided chemical equilibrium is closely 
approached. If the minerals are not binary, the concentration of other compo- 
nents because of the diadochic or substitutional relationships may affect the 
orderly distribution as discussed before. In fact, the approach to chemical 
equilibrium can be studied with respect to each component individually. Figure 
3.5 shows the distribution of Mn in coexisting minerals from charnockites 
(Saxena, 1968b). In other rocks as well as such orderly distribution of Mn is 
common. The distribution of Fe 2+ and Mg 2+ between coexisting olivine and 
orthopyroxene at 800°C and 900°C was experimentally studied by Medaris 
(1969). Although Medaris made repeated grinding and heating of the reaction 
products, Figure 3.6 shows that the distribution points both at 800 and 900 °C 
show some scatter. The difficulties are related to the kinetics of the ion-exchange 
reaction as equilibrium is approached particularly when the distribution approaches 
1:1 ratio in the two minerals. In contrast to these experimental results, the 
partitioning of Mg 2+ and Fe 2+ between orthopyroxene and Ca-pyroxene in meta- 
morphic rocks as studied by Kretz (1963) is remarkably orderly. Most distri- 
bution points fall on a smooth curve (Figure 3.7) and the distribution curve 
representing igneous rocks is clearly separated from the distribution curve for 
the metamorphic rocks. 
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Studies of partitioning of cations between coexisting minerals in natural 
rocks by several petrologists such as Albee (1965), Annersten (1968), Binns 
(1962), Butler (1969), Mueller (1960), Kretz (1959), Gorbatschev (1969), Hietnan 
(1971) among many others are attempts to rationalize the concept of metamorphic 
facies on a mineral-chemical basis. In experimental systems similar attempts 
have been made by Nafziger and Muan (1967), Larimer (1968) and Medaris (1969) 
among others. The results of such partitioning studies have generally confirmed 
the usefulness of the approach and the need for more thermodynamic data on 
crystalline solutions. 

In essence problems of phase quilibria are distribution problems and a 
statistical approach to such problems may be made to avoid the consideration 
of the thermodynamic properties of solution in individual minerals. Such ap- 
proaches have been made principally by Greenwood (1967) and Perry (1972) and 
should be applicable in solving the petrogenetic problem of incompatible as- 
semblages and the recognition of chemical equilibrium in natural or experimental 
systems. 


Table 3.1 

Composition of Coexisting Phases in Ternary Regular Solution 



a 





B 

B 


B 

A 

C 

B 

A 

C 

.010 

.082 

.908 

.043 

.868 

.089 

.108 

.047 

.406 

.040 

.079 

.881 

.174 

.718 

.109 

.336 

.195 

.478 

.080 

.071 

.849 

.339 

.522 

.140 

.530 

.394 

.576 

.100 

.066 

.834 

.413 

.431 

.156 

.602 

.489 

.632 

.130 

.056 

.814 

.512 

.309 

.179 

.699 

.624 

.715 

.170 

.037 

.793 

.625 

.171 

.205 

.821 

.785 

.796 

.200 

.020 

.780 

.700 

.081 

.219 

.90S 

.896 

.862 


r - B , _*g*A 
B B + A ’ 
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bution of Fe 2+ and Mg 2 + between synthetic olivine and pyroxene, 
bution data from Meolaris (1969). 




Section 4 


MEASUREMENT OF COMPONENT ACTIVITIES BY 
ANALYSIS OF PHASE DIAGRAMS 

An experimental measurement of activities of components in a crystalline 
solution particularly the silicates is beset with large difficulties, and the meas- 
ured values are subject to large errors. Therefore, it is an attractive proposi- 
tion to obtain such activity-composition relations from phase diagrams. 

However, here again there is no direct method of doing this without involving 
some kind of a solution model. The use of a solution model brings in certain 
uncertainties in the activity values which depend in extent and form on the 
choice of the model. The attempts to obtain the thermodynamic functions of 
mixing through the use of various solution models is still useful. For some 
crystalline solutions it may be possible to experimentally determine such 
properties. A comparison between the experimental values and the values 
based on a particular model would provide a greater understanding of the 
interrelationship of the crystal structural parameters on which the model is 
based and the thermodynamics of the crystal phase. In other cases where 
experimental determinations cannot be made, the empirically derived functions 
of mixing may be tested for their physical significance by their success in ap- 
plication to problems of petrogenesis. 

Symmetrical Mixtures 

Simple mixture 

The thermodynamics of a simple mixture or regular symmetric solution 
has been discussed before. Let us consider that the solution (A, B) M unmixes 
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into two coexisting solutions a rich in AM and / 3 rich in BM. At equilibrium we 
have 

3nd 

so that according to ’simple mixture’ model 

uf + R T in + W (1 - x%) 2 

= fiP° +R Tin x% + W ( 1 -x{) 2 
where <x Q and fi Q stand for the same pure endmember structure AM. Eliminating 

and/^° and substituting x& = 1 - * A a we have 

RT In + r (1 - x^) 2 = R T In (1 - *J) + IP (*“) 2 (4 * 3 ) 

or 


(4.1) 


(4.2) 


JL 

R T 


In {(1 - *?)/*?> 


(4.4) 


l-2xl 


This expression is similar to the one obtained by Thompson (1967). The equa- 
tion for the curve of coexistence of two phases may also be written in terms of 
critical temperature T c of unmixing and the mole fractions by substituting 

W - 2 RT- 


in (4.4) 


T = 2 t A (4.5) 

C .ln {(1 -*“)/.£ 

If we have data on the composition of coexisting phases at different temper- 
atures and the form of the solvus is symmetric, the value of w and the activity- 
composition relations can be calculated. 
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Symmetrical mixture of higher order 


Symmetrical crystalline solutions may not be simple mixtures and require 
a two or more constant expression for representing G m as follows: 

<4,6) 


For such a solution the relation (4.3) is 


R T In *J + {^ + Aj (1 - 2 *“) 2 } (1 - *“) 2 

= R T In + (A 0 + A 2 (1 - 2 x.J) 2 } (1 - xP) 2 


(4.7) 


Using the relation^ “ ={j*P , we can write a similar relation as (4,7) and then 
solve the two equations simultaneously to obtain A Q and A 2 and the activity- 
composition relation 


Asymmetrical Solutions 

Sub-regular model 

As mentioned before the excess free energy of mixing may be expressed as 
a polynomial in the mole fraction * A or * B for the compound (A, B) M according 
to Guggenheim's equation: 

^EM = *A X B ^0 + 4 1 ( X P> ~ X b) + 4 2 ( X k ~ *b) 2 + * ^ ^ 

If we put A 2 - 0 in (1.47) we have a two constant equation for an asymmetrical 
solution. Proceeding as in the previous sections we have 

R T In + R T In /J = RT In + R T In fP (4.8) 

R T ln*“ + R T In = R T In x? + R T In < 4 - 9 ) 

B D D D 

Substituting values of RT Inf from (1.48) and (1.49) we have 

RT In X* + (x£) 2 {A Q + A t (3 x\ - xj) = R T In + (xg) 2 {A Q + (3*J-*g)}( 4 *l<>) 
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(4.11) 


RTinx J + (*“) 2 ty, + A, (3 - *;» = R T In + (xft* 

(A Q + A x (3 4' x^» 

The two independent relations (4*10) aind (4.11) now can be solved simultane- 
ously to yield the values of the two constants A and A t . 

The method of calculation presented above is equivalent to that used by 
Thompson (1967) and Thompson and Waldbaum (1969), Thompson's (1967) 
equation for the excess free energy of mixing is 

C - x G + x„ G (4.12) 

°EM X A u 2 + *B U 1 


where 


and 


G„ = 


*A W G2 


G, = 


X A *B *G1 


This is as if the crystalline solution is composed of * A moles of a simple mixture 
with W G2 and * B moles of another simple mixture withlP G1 . We have then 

^EM _ (*G1 *B + *G2 X a) X A X B 

= »gi (1 - **) + F 0 J u - *b» *a *b 


„ ( 2 - 2 *a) „ <2-2» b ) 

*gi j +, ° 2 2 *a*b 


(4.13) 


substituting 1 = ^ in (4.13), we have 


W. 


g em “ 


G1 


V. 


^d-*A + *b)+4 2 (1+«A-V- 


X A X B 




*G2 + *G1 *G2 “ *G1 , N 

9 + 9 ( *A “ *B> 


X A X B 


(4.14) 
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which is of the same form as Guggenheim's equation with two constants A and 
A V 


Therefore 


9' + W W C7 - If' 

/L = l 1 and A. — 

0 2 1 2 


(4.15) 


A 0 /R T and A x /R T would correspond to the notations B Q and C G respectively 
used by Thompson (1967) following Scatchard and Hamer (1935). 

For the activity coefficient, we have 

R T ln f A = O b ) 2 lA 0 + A 1 (3 * A - x B ) + . . .] (1-48) 

Substitution of relations (4.15) in (1.48) gives 




*02 + *01 *02 "*01 


(*A - X B + 2 *A> 


- *-2 


= x 


B 


W 4 - w w - w 

n G2 ■ *G1 h G2 H G1 

+ 


(*A - X bH + ( ff C. 2 “ *fll> * 


g2 G1 t A 


= x 


B ^ ^G1 X B + ^G2 ^A^ + ^G2 “ W CA *a)^ 


G2 A G1 A' 


- ~2 


X B ^G1 (*B “ *A> + ^ ^G2 
_ X B ^G1 (1 " 2 * A ) + 2 W Q2 * A ] 


X B ^G1 + ^ *A ( ff G2 “ ^Gl^ 


which is the same as used by Thompson (1967). 


(4.16) 


Substituting the calculated values of A 0 and A 1 in (1.47), the excess free 
energy of mixing can be estimated. We can then determine the solvus bounding 
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the two phase region by the graphical tangent inethod or by a suitable iteration 
numerical method. The calculated values of “ , x@ and are then compared 

ABA B 

to the observed mole fractions to test the applicability of the model. Other 
excess functions of mixing may be calculated by using the relations (1.50) and 
(1.51). 


Quasi-ch emical approximation 

This model has been discussed before. Green (1970) used the model to 
study the halite-sylvite solvus. Let us consider two coexisting phases M and 
N with components A and B. We have at equilibrium at a certain p and / . 

/Li” -fj . N and 

or 

[j. m + R T in xf + ft T In /jj = /x AN + R T ln xj + R T In f* 
and 

/x BM +RT In xg +RT In fM + RT In xN -|. RT in fg 

Since both M and N obey the same equation of state the chemical potentials of 
pure A in M and A in N and B in M and B in N are cancelled. 

Substituting values of / and f from (1.33) and (1.34) in (4.17) we have 

A 6 


In ^ + — — 1 n 1 1 + 


0g 00- T) 


05 0 + i)j 


I z q. 

= In*^ + ln 

A 2 


l t 


*5 0' - 1)" 


(4.18) 


L $1 0 ' + i)J 


Z O, 

ln + - In 

B 0 


1 + 


- 1 )' 

4>l (fi' - 1) 


z 7, 

In f 1 In 

B 2 


1 + 


^ 0' - 1)“ 


(4.19) 


L B 


tfg c/ 3 ' + l) 
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where P and/3' correspond to phases M and N respectively and along 0’s are 
defined by equations (1.28 using instead of % k etc.) and (1.35) presented be- 
fore. q andcy are contact factors discussed before. They are not independent 
and should approach unity simultaneously. Green (1970) assumed/^ = 1. 

The two independent relations (4.18) and (4.19) contain two unknowns q /q o and 
W and can be solved by an iteration process. The ratio q l /q 2 is a function of 
the geometry of the substituting chemical species and therefore may be regarded 
as almost independent of temperature. Substitution of q /q back in equations 
(4.18) and (4.19) gives two independent values of W at each temperature. Any 
difference noted in the two values of W would be due to the inadequacy of the 
solution model to fit to the experimental data. 

The solvus bounding the two phase region may as before be determined 
graphically by the double tangent method on a plot of free energy of mixing 
against mole fraction or by a numerical iteration method. The excess functions 
of mixing can be calculated using (1.36) and (1.37). 

Calculation of Functions of Mixing Example: 

The System CaW0 4 - SrWC> 4 

Chang (1967), presented the data on the two phase regions with solvus in 
the binary tungstate R 11 W0 4 type crystalline solutions. Table 4.1 shows 
the data on the composition of the coexisting phases a and/3 rich in CaW0 4 and 
SrW0 4 respectively. Calculated values of A Q /RT and A 1 /RT according to the 
subregular model are listed in Table 4.2, The relationship between A Q /A 1 and 
T is linear and given by 

A 0 = 22212.0 - 16.7278 T < 4 * 20 ) 
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(4.21) 


A 1 = 3476.67 - 4.4379 T 

We may then calculate excess free energy of mixing from the relation 

^EM " X SrW0 4 *CaW0 4 ^0 + A 1 ( X SrW0 4 “ X CaW0 4 ^ 

and the activity coefficients from the relations 

R T in f SrWo = A'caWO ^0 + ^1 ( 3 *SrWO " *ChV'0 
4 4 4 4 

^ T ln ^CaW0 4 = ^SrW0 4 ^0 " /} 1 ( 3 X CaW0 4 “ *SrW0 4 ^ 

Figure 4.1 shows the activity-composition relation at 800 and 1000°C. We 
may also plot free energy of mixing against composition and find the 'tangent 
method* the composition of the coexisting phases. In the present case the dif- 
ferences between compositions calculated by the model and the observed com- 
positions in Table 4.1 are found to be small. 

As mentioned before the activity-composition data and other thermodynamic 
functions as calculated from phase diagrams are sensitive to the nature of the 
assumptions and the model used. For the system CaW0 4 - SrWO , it is found 
by actual calculations that the use of quasi-chemical approximation predicts 
solvus with similar accuracy as does the sub-regular model. The use of quasi- 
chemical approximation requires the values of q , q 2 and the coordination z. 

z is the number of the nearest Ca or Sr ions surrounding each other. In CaWC> 4 

0 

there are four Ca 2 surrounding each Ca 2+ at distance of approximately 3. 9 A. 

o 

There are four more Ca 2+ at a distance of 5 A. We may assume 2 as 4 and solve 
the two equations of the quasi-chemical approximation (4.18 and 4.19) simultane- 
ously to find q /<? . 
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In the present case, however, we may assume either that q x + q 2 - 2 or 
A i q 2 = 1. The differences in the calculations of W using either of the two as- 
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sumptions are small (see also. Green, 1970). A computer program now may be 
used which solves each of the two equations independently by using various values 
for q 1 and q 2 and compares the IF values so obtained until it finds the best set of 
IF values which match each other. Such values with Z as 4 are listed in Table 
4.3. q x and q 2 are 1.20 and 0.80 corresponding to SrW0 2+ and CaWO 2 * 
respectively. The atomic radii (K) for Sr 2+ and Ca 2+ are 1.12 and 0.99 (Ahrens, 
1952). We have then q 1 /q 2 as 1.50 and B (Sr 2+ )/K (Ca 2+ ) as 1.13 which are 
not similar. 


The following equation describes the relation between the calculated if 
( average of the two values listed in Table 4.3, second column) and T : 


R T- 34.64 - 0.05506 T + 0.000025 T 2 


(4.22) 


Figure 4.2 shows a comparison of G m at 800 °C calculated according to the 
sub-regular model and according to the quasi-chemical model. The excess 
free energy of mixing according to the latter is nearly twice of that calculated 
according to the former. Differences between the other calculated functions of 
mixing// and s n . are even more marked (Figure 4.2). Unfortunately we do 
not have data on experimentally determined //_„ and S_, w for the CaWC) -Sr WO, 

-KM IlM 4 4 

system and, therefore, no way to know which model predicts the thermodynamic 
functions better in this particular case. 


For the system NaCl-KCl, Green (1970) provided the comparison in the 
thermodynamic quantities calculated by using the sub-regular model and the 



■-I 


■< 


'•& ■ 
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quasi-chemical model and those determined by experiments. The thermodynamic 
quantities predicted by quasi-chemical model are closer to those measured ex- 
perimentally. 

A comparison of the predictions of the functions of mixing in several binary 
alloys by regular solution model and by the quasi-chemical model (Lupis and 
Elliott, 1967) shows that generally the predictions by the latter for the excess 
free energy are closer in agreement with experimental determinations than 
those by the former. The prediction for the excess entropy by the quasi-chemical 
model are not satisfactory. This may be in part due to the neglect of the non- 
configurational excess entropy in many of the binary alloys. For the halite- 
sylvite system, Green (1970) finds that the nonconfigurational contributions are 
unimportant and suggests that the positive excess entropy of mixing found in the 
NaCl-KCl system may result from the introduction of vacancies or other defects 
into a crystalline solution. 

This approach of calculating thermodynamic functions of mixing by the 
analysis of phase diagrams is relatively new in. the field of mineralogy and 
deserves more attention from the mineralogists and petrologists. The fact 
that there is no unique analysis of a solvus and two ore more solution models 
may be applicable to the same solvus data need not deter us from acquiring 
more phase data and interpreting the same with the help of various solution 
models. Experimental verification of many of these results may not be possible 
in the near future. However it may be possible to test such thermodynamic 
data by their application to petrogenetic problems and by getting in return some 
physically meaningful results. 
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Coexisting Phases with Different Crystal Structures 
In minerals such as pyroxenes and amphiboles , one often finds a pair of 
minerals whose composition lies on a binary compositional join but whose 
structures differ more or less from each other. A case in point is the enstatite- 
diopside pair. The compositions in this series lie on the binary MgSiO -CaSiO 

o J 

join. The solvus on this join is not comparable to the solvus for the NaCl-KCl 
system or for the CaW0 4 -SrW0 4 system. The important difference is that 
enstatite and diopside do not obey the same equations of state as do the NaCl 
and KC1 or CaW0 4 and SrW0 4 . Phase equilibrium data on the enstatite-diopside 
solvus such as that of Boyd and Sehairer (1964), therefore, cannot be used for 
evaluating the activity-composition relations by the methods of previous sections. 


Let us consider two coexisting phases a and j3 with different structures but 
the same composition (A, B) M. We have 


< 4 ' 23 > 


and 


LL a = U. P 

^bm m 


(4.24) 


where 




a \ 2 


(4.25) 


and similarly others. 


Substituting relations such as (4.25) in (4.23) and (4.24), we have 


A + RT In + W* (1 - x^)2 = + rt In (1 - 


(4.26) 


M BMa + R T In + w° (1 - X^)2 fe + R T In + W^ 3 (1 - x^)2 


(4.27) 
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Since a and /3 do not have the same structure and, therefore, do not obey the 
same equation of state generally ^ AMa 4 and p. BMa As such we need 

to know the values of these chemical potentials at the given P and T to be able 
to solve the two equations (4.26) and (4.27) simultaneously to find and WP , 

The free energy diagram in Figure 4.3 shows this situation. In a binary system 
at a given P and T % we may have two coexisting phases with compositions repre- 
senting the minimum of free energy. 

The enstatite-diopside system may be similar to the example shown in 
Figure 4.3, The binary join is (Mg, Mg) Si 2 O g - (Ca, Mg) Si 2 O g . There are 
further complications here since (Ca, Mg) Si O with orthopyroxene structure 
and (Mg, Mg Si 2 O g ) with diopside structure are unknown. 
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Table ’4.1 


Chemical Composition oi Unmixed Phases in the System CaW0 4 -SrW0 4 


T°C 

'Sr 


*5r 


550 

0.005 

0.995 

0.995 

0.005 

600 

0,010 

0.990 

0.980 

0.020 

650 

0.025 

0.975 

0.955 

0.045 

700 

0.035 

0.965 

0.905 

0.095 

750 

0.067 

0.933 

0.800 

0.200 

800 

0.120 

0.880 

0.630 

0.370 


a and /3 are coexixting phases rich in CaW0 4 and SrW0 4 respectively. The compositions are 
from Figure 3 in Chang (1967). 





Table 4.2 


The Calculated A n /RT and A Jrt In (Ca, Sr) WO 

U 1 4 


T°C 

A 0 /RT 

A t /RT 

600 

4.411 

-.278 

650 

3.596 

-.269 

700 

3.113 

-.455 

750 

2.537 

-.514 

800 

2.015 

-.611 

900* 

1.106 

-.743 

1000* 

.352 

-.860 

* From equations (4.20) and (4.21) 
Note: An error of d5 percent in the mo 

e fractions (Table 4.1) results in 

d200’cal/mole error 


in determining A Q and Ay 



Table 4.3 

For the System CaWO . -Sr WO, 

z R T 4 4 















Section 5 


MEASUREMENT OF COMPONENT ACTIVITY USING 
COMPOSITION OF COEXISTING MINERALS 

Experimental data on the distribution of a component between two coexisting 
crystalline solutions at a fixed p and T for systems such as olivine and pyroxene 
have been collected by NafzXger and Muan (1967), Larimer (1968) and Medaris 
(1969). Distribution data are also available for natural assemblages but the P 
and T of their formation are indefinite. The distribution data from natural 
assemblages in many cases may be found to represent ion-exchange equilibrium 
closely. If we are not concerned with precise P and T, we may use such data 
to obtain useful information on the thermodynamic nature of mixing in the 
minerals. For this purpose, the thermodynamic equations according to various 
solution models for binary solutions presented in this section may be used. 

The composition of coexisting phases which do not obey the same equation 
of states may be used to find the activity- comDoeition relations in each phase 
in suitable cases. Consider a and /3 with chemical formulae (A, B) M and (A, B) N 
respectively which are in ion-exchange equilibrium at a certain P and T , We 
have 

Aa+B/3^Ba + A/3 ( 5a ) 

The equilibrium constant K 5tt is given by 


The term in the first bracket is the distribution coefficient K^. Depending on the 
nature of the data available, the following cases may be considered. 

Compositional Data Available on a Complete Distribution Isotherm 
We have the choice among simple mixture model, the two constant asymmetric 
model and the regular solution model with quasi-chemical approximation. As 
discussed before, according to simple mixture model, we have 


W a 


ln *,. = ln*D-jf_<l-2*J) 


+ — (1 - 2 x%) 

RT k a ' 


(3.8) 


A non-linear least squares fit using the data on x* and xP finally yields K. , 

A A 5 8 

W a and W?, 

According to Redlich and Kister equations (King, 1969, p. 326) we have 


R ‘ ln = *B H, + *1 < 3 *A - *B> + Aj (*A - *„) <5 * A - V + >, < 5 ' 2 > 


R T ln f s = x l {A o - A 1 (3 * B - * A ) + A, (*„ - , A ) (S * B - « A ) +-. .} 


(5.3) 


and therefore 


R T ln * = A a (*„ - * A ) + 4, (6 * A Xg - 1) + A 2 (» B - » A ) (1 - 8 * A * B ) 

B 


(5.4) 


Substituting the values of f£/f? and fP/fP by using (5.4) in (5.1), neglecting 

O ^ A *5 

the constants A 2 's and rearranging in logarithmic form we have 


In K = 


A a 

n o 




*- sln * D + jir(‘ I ‘*5 )+ jiT (6 "b^a- 1 ) 






+ — — (*£ - *5) + — (6 x^ x£ - 1) 
r T K B ^ R T K Ab ' 


(5.5) 


The equation (5.5) is of the form 


N = M + Aj Xj + A 2 x 2 + A 3 x 3 + . 

where N and x’s are known quantities. It may be solved by a numeric least 
squares method yielding M(= In X ga ) and other constants. We must have a mini- 
mum of five distribution points. 

According to the quasi-chemical approximation, f A is given by 



(P “ 1) 

_ 0 + 1) 


zq A /2 


(1.33) 


and similarly L as in (1.34). Substituting these values of f *s in (5.1) we have 


1 + 


. *NZ a tr a /2 


4>i (P* + 1>. 


1 + 


4£ 


V r B v ■ 

V ~ 

C , -K ‘ - 

1 ^(p“+i)J 


^ (Z^ + 1). 


1 + 


K - 1 ) * 

^ (Z 3 ^ + 1). 


B B„ r = ^Sa 
z H q%/2 


(5.6) 


where Z 3 and Z 3 ' are for phases a and /J respectively and are given by (1.28) with 
x A replaced by 0 A etc. and ^ are contact factors for phases a and /3 
respectively. A numerical least squares method may be used to solve (5.6). 
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Compositional Data Available on a Complete Distribution 
Isotherm And The Activity-composition Relation in 
One of The Two Coexisting Phases 

Depending on the accuracy of the data, we may use equations (3.8), (5.5) 
(5.6). If necessary we may use all the three constants in equation (5.4). We 
may write (5.4) as 


RTln-jr- =X B ( 2A 0 + 6A l +1 ° A 2> ■ X B ( 6 A 1 + 24 A 2> +X B ( 16 A 2^( A 0 +A 1 + A 2^ 

' B // 


(5.7) 


Transforming (5.1) in logarithmic form and substituting values for fP /fP 

A B 

from (5.7) we have 


A n + A, + A, 


InK^ + lnH = (lnK s , + “° T '‘ lT " 2 ' ' ‘ 


n 


R T 


- x 


2 An + 6 Aj + 10 A 2 


B' 


R T 


+ X 


'6 A t + 24 A\ 


B 


R T 



(5.8) 


Equation (5.8) is of the form 

Y = a Q + a x x + a 2 x 2 + a 3 x 3 + . . . 

If we know the activity-composition relation in a, we may solve (5.8) by least 
squares analysis. 

If sufficient number of distribution points are available one may use (5.5) with 
the third constant A 2 and compare the results with those obtained by using (5.8). 
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As an example of calculating the activity-composition relation in binary 
solutions using the distribution data, we consider the olivine-chloride solution 
system. 


Activity-composition And G m in Olivine (Mg 2 Si0 2 — Fe 2 Si0 4 ) 

Schulien, Friedrichsen and Hellner (1970) studied the distribution of Fe 2+ 
and Mg 2+ between a 'chloride 1 solution and olivine at 450 ° to 650 °C. These 
data may be used for calculating the activity-composition in both olivine and 
in the 'chloride* solution. The nature of the 'chloride' solution in the experi- 
ments may not be clearly known. This need not concern us for the present and 
we assume that it is a homogeneous solution phase and present the ion-exchange 
equilibrium as done by Schulien et. al. as 

FeCl 2 +iMg 2 Si0 4 ^ MgCl 2 + |Fe 2 SiO, (6b) 

solution olivine solution olivine 

The equilibrium constant is given by 


x ii fr 1 

^ _ Mg Fe Mg Fe 

5b v sol v ol fsol col 
Mg l Fe r Mg 

We may now use one of tlje three methods represented by equations (3.8), (5.5) 
and (5.6). In the present case, it was found by actual computations that consistent 
thermodynamic parameters for the isotherms at 500, 550, and 650 °C are given 
by only the 'simple mixture' model. Obviously the date, accurate as they are from 
an experimental point of view, are not accurate enough for needing a two or more 
constant equation (5.5) or the quasi-chemical approximation (5.6). 


W /RT values calculated by using the Simple mixture' model are listed in 
Table 5.1. The data at 600 °C was not considered because of several scattered 
points. It is noted that W /R T for olivine at 450° C is greater than 2.0 and is 
unrealistic. If we consider a linear relation between ff /o1 and l/T (Figs. 5.1 and 
5.2), we find that W ol /RT at 450°C is 2.0 and similarly W sol /R T and K sb are 
0.91 and 0.76 respectively. 

Figure 5.3 shows the distribution isotherms calculated by using the data 
in Table 5.1 at 500, 550 and 650°C. The isotherm at 450°C is plotted using the 
linearly extrapolated values of K sb , [f /o1 and ff /so1 . The model does not fitto these 
distribution data for the Mg-rich chloride solution and olivine. Figure 5.4 shows 
the activity-composition relation in olivine (Fe 2 Si0 2 — Mg 2 Si0 2 ) at 500 and 550° C. 
At 650 °C olivine is close to ideal solution. 
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Table 5.1 

Inter-change Energy IT and the Standard Free Energy Change AG° 
in Olivine-chloride Solution System 


T°C 

IJ/sol 

FT 

W oL 
R T 

-RT/nK 

= AG° . cal/mole 

5b ' 

450 

1.86 

2.32 

1433 

500 

0.87 

1.28 

723 

550 

0.59 

0.73 

421 

650 

0.52 

-0.014 

181 
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Figure 5-4. 'Activity-composition relation in olivine 
At 650 °C, olivine is found to be nearly 
Schulien et al. (1970). 
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Section 6 


ORDER-DISORDER IN Fe 2+ -Mg 2+ SILICATES 
Long-range order-disorder phenomena in ferromagnesian silicates differs 
from that in alloys in several important respects. First, as opposed to alloys, 
only a certain number of cations take part in the site exchange while the silicate 
framework remains more or less inert. Second usually a complete crystalline 
solution exists between the Mg and Fe end-members. Third since Fe 2r and 
Mg 2+ are similar in size, charge and other characteristics, the site preference 
energies (corresponding to the difference in binding energy of the ion between 
the nonequivalent sites) arc not strongly dependent on the degree of order as is 
usual in many binary alloys. Order-disorder or the intracrystalline cation 
distribution in silicates is measurable by x-ray (see Ghose, 1961) and other 
spectroscopic techniques. The energy of the intracrystalline ion-exchange is 
part of the Gibbs free energy of the crystal and is, therefore, a very useful 
thermodynamic quantity (see Mueller, 1969, Thompson, 1969). 




where % refers to the mole fractions , K a and K „ are rate constants and are 

pa ap 

fimctions only of P and T , 4>s are analogous to activity coefficient products in 
macroscopic chemical system and are functions of P , T and</>. At equilibrium 
we have 


and 


d *A 

dt 


= 0 


~\ 


l 

V 



K 


6a 


_ V „ *1 fl f* 


K 


a/3 

^ B ■' B A 1 A 


a B a l 


( 6 . 2 ) 


where f is the 'partial* activity coefficient and a the 'partial' activity. The 
product of the f s appear as 0in (6.1). The term 'partial 1 is used to distinguish 
between the activity of A on the site from the activity of A in the crystal. 

The distribution coefficient is 


*D = 


* A 

A B 


a 


*B*A 


& D has sometimes been referred to as the ordering parameter. The distribution 
coefficient, however, should not be confused with the ordering parameter S used 
to describe ordering in alloys. S equal to unity corresponds the highest possible 
order and S equal to zero to complete disorder. This is opposite in the case of 
K^. Further k b will be used to describe order-disorder in non-stoichiometric 
silicates forming complete crystalline solution series. In such silicates the 
formation of a fully ordered or disordered periodic structure is not possible. 
Even with the greatest tendency towards ordering some of the excess atoms of 
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one component must inevitably occupy sites belonging to the other which leads 
to a lower order or disorder. The distribution coefficient is a function of T and 
the varying A:B ratio in the crystal and, therefore, it is of little thermodynamic 
significance. 

The equilibrium constant K- 6a is a function of P and / only. However, as the 
volume changes involved in the ion-exchange are negHgible, we ignore the de- 
pendence of K on P and consider it only as temperature dependent. 

The definition of chemical potential of a cation on a site presents certain 
problems (Mueller, Ghose and Saxena, 1970). One may write, as done by Grover 
and Orville (1969), for the chemical potential of a cation A on the sites a and /3: 

y-l = + * T In al (6 ' 3) 

1*1=^+ R T lnog < 6 - 4 ) 

where the u® -s are the standard chemical potentials, and the a ! s are the cor- 
responding ’partial' activities. According to classical thermodynamics, how- 
ever, it is incongruent to define two different chemical potentials for one species 
in a single homogeneous phase. In such a case we have: 

and =^a° 

To avoid the above difficulty Borghese (1967) regards A in site a as a dis- 
tinct species from A in site/3. This is somewhat analogous to speaking of the 
chemical potentials of 0 2 and 0 3 in a homogeneous gas phase. The idea of de- 
fining a new potential analogous to chemical potential called a site preference 
potential (Greenwood in Grover and Orville, 1970) could also be considered. 


The use of a 'site preference potential* may be avoided in practice, partic- 
ularly as we are rarely concerned with its quantitative use in thermodynamics. 


The standard site preference energy or the intracrystalline ion-exchange 
energy AG° for the reaction (6a) is given by 

A G° = - R T In K a < 6,6 ) 

where K a is the equilibrium constant and is a function of T only unlike which 
is a function of both T and composition. 

Thermodynamic Functions of Mixing 

One of the principal aims of the study of order-disorder phenomenon is to 
investigate the thermodynamic properties of the crystalline solution as a whole. 
In case of an ideal macro-phase the activity- composition relation is given by 

(6.7) 

a i = ( x i) N 

where N is the number of structural sites in the crystal. In the case when 
there are two sites we may have 

“i = (*i) <*i> 

and if these sites are different we can write 

°i = {(*i + *f)/2> 
or 

°i =<*i) a (*i)^ 

where a and ft are two nonequivalent structural sites, 
been generally used (Mueller, 1962, Thompson, 1969). 
method to the non-ideal case, we may write 
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( 6 . 8 ) 

(6.9) 

The latter method has 
Extending the above 



( 6 . 10 ) 


«1 = K) a (“i) /3 

where a ® andaf are 'partial' activities referring to the sites. In an orthopy- 
roxene (Mg Mg Si 2 0 6 - FeFeSi 2 0 6 ), where we have two sites Ml and M2, the 
activity of Fe 2+ in the crystal may be expressed as 


a °P X — o M1 a M2 
Fe 2 + pe 2+ Fe 2+ 


( 6 . 11 ) 


If we consider the activity on one cation basis i. e. for the crystal (Mg Si0 3 ) 
Fe Si0 3 ) we have 


a ° P x _ (o Ml_ 


M2 


Fe 


2 + 


Fe 2+ Fe 2+ 


y/2 


( 6 . 12 ) 


The 'partial' activity af is equal to f L a where/ is the partial activity coef- 
ficient. At a certain temperature the atomic ratio*, in the two sites a and f3 
can be determined by x-ray or other resonance techniques. The next problem 
involves the evaluation of the 'partial' activity coefficients. 


Several crystals of suitable composition (A, B) M between the endmembers 
AM and BM may be chosen and heated at a certain temperature for time long 
enough to attain equilibrium for the intra-crystalline ion exchange (6a). Several 
such distribution isotherms may be obtained. We shall be required to find a 
suitable model for interrelating the partial activity coefficient with the atomic 
fraction at the site. The ’simple mixture' or regular solution model may be 
found useful in cases where the form of the distribution isotherms do not indi- 
cate two much non-ideality of mixing A and B at a.and/3 . We may write 

In K a = In K D - (1 -2*J) + Ii(l -2*£) (6.13) 
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where IP is related to the partial activity coefficient by 


Rfln/J=r(l-*;) a 


(1.25) 


At this point certain other partial functions of mixing may be considered. 
The partial free energy of mixing at the sites is given by 

(6.14) 


GZ = xlRTlnal + x*RTlna* 


= R T In fl + R T In f* + T (*J R In R in *£) 


(6.15) 


Substituting S * s -R (*“ ln * * + x* In *“) and (1.25) in (6.15) we have 

IV1 A n D O 


~ W * *A "B-^M 


(6.16) 


The term W a xf . is also the partial excess free energy of mixing. 

A B 

In a crystalline solution such as AAM - BBM the total free energy of mixing 
is given by 


G m R T ln a Z"< + (- S l^) S T ln “b'°b 


(6.17) 


which can be finally shown to be 


g m + C M 


(6.18) 


„a. „ B 

x . — x\ 


RT ln K 6b - T (5® + Sg) + JP a ' 6 * 19 > 


This is a similar expression as derived by Grover and Orville (1969) for ideal 
mixing at the sites. Note that the above expression £ 6a is equilibrium constant 
and not the distribution coefficient as in the case of ideal mixing. 


72 


Substituting A G°= - R T lntf in (6.19) we have 

G m=^T^ AC °- 7 '( S S +S S) + ': C ™+G^> <6.20) 


Thus the free energy of mixing in the crystal as a whole is a result of a) 
energy due to the distribution of the cation between a and /3 sites b) entropy 
change due to the distribution of A and B within a and /3 sites and finally c) the 
excess energy of mixing arising due to the non- ideality of solution in a and/3 , 


The partial excess free energy of mixing at the sites are 


pa. _ ii a 'poo- 
^EM “ ,? EM “ 1 0 EM 


( 6 . 21 ) 


and 


gP = hP - T sP 

EM n EM 1 "EM 


( 6 . 22 ) 


Substituting (6.21) and (6.22) in (6.20) we have 




M 


A G» + (HJ, ■ + Hi) - T <8? + + Bf + Sfc) 


(6.23) 


The use of the thermodynamic relations presented in this section and before 
will be made in analyzing the data on site occupancies in orthopyroxene later. 

Kinetics of Order-Disorder 

Virgo and Hafner (1969) made an important observation that there is an 
apparent cut-off or transition region on the temperature scale below which no 
more ordering or disordering occurs.. This transition temperature in orthopy- 
roxene was estimated to be approximately 480° C. Above this temperature the 


activation energy required for diffusion to start in the direction of disordering 
is of the order of 20 kilocalories (Virgo and Hafner, 1969). Below this tempera- 
ture the activation energy should be very high. This is confirmed from the 
measurement of order-disorder in metamorphic pyroxenes which cooled slowly 
through geological time. Figure 6.1 shows the data on the K D values for the 
distribution of Fe 2 + and Mg 2+ between Ml and M2 sites in metamorphic ortho- 
pyroxene. From these data and the experimental results shown in Figure 8.4 
later, it is noted that no orthopyroxene shows a degree of order representing 
temperatures lower than 450° C. 

Mueller (1970) proposed a two-step mechanism for order-disorder kinetic 
in silicates. This involves a low temperature process with high activation 
energy and a high temperature process with a lower activation energy. This 
mechanism may be responsible for ordering characteristics distinguishing 
metamorphic, igneous plutonie and volcanic pyroxenes. The intracrystaiiine 
ion-exchange equilibria in igneous plutonie rocks is not ordinarily quenched in 
at any temperature because of slow cooling rate. The same applies to such 
equilibria in metamorphic rocks. However in these rocks attainment of such 
equilibria is possible below the transition temperature if crystallization or re- 
crystallization occurs at these temperatures. In volcanic rocks because of 
rapid cooling, there is a high possibility that the intracrystalline equilibria is 
quenched in and the temperature indicated by order-disorder is not very much 
lower than the original temperature of crystallization. 


500 


• « • 


POTENTIAL BARRIER 


- 1.5 


Figure 6-1. Temperature estimate of the ion-exchange equilibrium 
between sites in natural orthopyroxenes. 450 °C ap- 
pears to be the transition temperature below which no 
more ordering takes place due to a'potential barrier. 
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ORTHOPYROXENE CRYSTALLINE SOLUTION 

Orthopyroxene is one of the few important rook forming minerals which can 
be considered as quasi-binary without significant loss of accuracy. Usually 
more than 95 percent of the mineral is a crystalline solution of the endmembers 
enstatite (MgSi0 3 ) and ferrosilite (FeSi0 3 ). Fe 2 * and Mg 2+ are distributed 
between two nonequivalent sites Ml and M2. With the use of X-ray or Mossbauer 
spectroscopic technique, it is possible to determine the proportion of Fe 2+ in the 
two nonequivalent sites (Evans, Ghose and Hafner, 1967). These data can be 
used with the help of suitable solution models to determine the thermodynamic 
properties of the solution (Virgo and Hafner, 1969, Saxena and Ghose, 1971). 

Inter-site Ion-exchange 

Thermodynamics of the intracrystalline ion-exchange was discussed in the 
preceding section (6). Ml and M2 sites may be regarded as two interpenetrating 
subsystems, each with its own thermodynamic properties of mixing. In analogy 
with heterogeneous ion-exchange equilibria, we may write the ion-exchange 
reaction as 


Fe 2+ (M2) + Mg 2+ (Ml)=^Fe 2+ (Ml) + Mg 2+ (M2) (7a) 


The equilibrium constant for the above reaction at a certain P and T is 


K , , = 

7a 


V M1 fMl vM2 fM2 
*Fe Fe *Mg *Mg 

V M2 CM2 v Ml fMl 
iTe A Fe *Mg *Mg 


(7.1) 
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K 7b is mainly a function of T . P has little influence. The standard free 
energy for ior-exchange at a certain 7 is 

A G 7a - - R T In K la (7.2) 

we have noted before that this energy is part of the total Gibbs free energy of 
the crystal and is, therefore, an important thermodynamic quantity. To determine 
AC 0 or K , we need first, the determination of the atomic fractions * M1 and 
and second, the determination of the partial activity coefficients f’s which 
are functions of T and composition. While the atomic fractions or the site 
occupancies can be determined quantitatively by using X-ray or Mossbauer tech- 
nique, the partial activity coefficients cannot be determined without the use of a 
certain solution model. 

Order-disorder on Individual Sites And The Choice of a Solution Model 
Figure 7.1 (Figs. 1 and 2, Ghose, 1965) shows a scheme of ordering in 
orthopyroxone crystal as a whole. In hypersthene (Fs so En so ) we have the possi- 
bility of complete occupation of Ml sites by Mg 2+ and of M2 sites by Fe 2+ . For 
all other compositions, one of the two ions must inevitably occupy the sites other 
than they normally prefer. In fact, as suggested before we except that the two 
ions will always show some kind of equilibrium distribution over the two sites as 
a function of temperature and composition. Although hypersthene has the right 
composition to be completely ordered, kinetics of the ion-exchange below a certain 
transition temperature precludes the ordered structure shown in Figure 7.1. 

For considering the ion-exchange (7a) we are required to consider Ml and 
M2 individually as sub-systems. Complete order in the crystal as a whole as 
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shown in Figure 7.1 also means complete order on the sites themselves. For all 
other compositions, we may consider whether two neighboring Ml sites or M2 
sites are both occupied by the same cations or by two different cations. We may 
assume that the occupancy of Ml and M2 sites by Mg 2+ and Fry*' is disordered 
and that the solutions at the two sites approximate to ’simple mixture’ model. 

We shall later consider the effect of ordering on the sites by using the quasi- 
chemical approximation and compare the results. 

Determination of Site Occupancy in Heated Orthopyroxenes 
By using MOssbauer spectroscopy, Virgo and Hafner (1969) determined a 
distribution isotherm at 1000 °C. Saxena and Ghose (1971) similarly determined 
isotherms between 500 and 800°C (Fig. 7.2). Orthopyroxene crystals are heated 
for a length of time sufficient to bring about the ion-exchange equilibrium. The 
samples are then quenched and the site occupancy is determined. Note that the 
site occupancy is averaged over all the sites present in the crystals and v, a 
function of the changes that take place not only in the neighboring cation but 
also in the remaining Si-O structural framework. 


Sites as Simple Mixtures 
Equilibrium Constant, IF M1 and W K2 

The distribution isotherms can be represented by the following equation 


wMl 

In K,. = In K, t' (1 - 2 *“) (1 - 2 *“) 


where K n is *“* (1 - x“;J)/x^ ( j _ „>»), 




* Fe the site occupancy and if’ the energy constant of the simple mixture model. 

A non-linear least squares program is used to calculate # 7a , IP Mi and IT M2 
by simultaneous iterations on all the three unknown variables. Figure 7.3 shows 
a plot of A C° Q the standard free energy of ion-exchange (= - RT In ^ ), lf ' M1 
and ir M2 against l/T . There is a linear relation between l/T and all the three 
variables at 600, 700 and 800°C. The equations to the three straight lines are: 


W M1 = 3525 f— j - 1667 


W M2 --- 2458 - 1261 


AC® = 4479 - 1948 (— 

a V J 


(7.3) 


(7.4) 

(7.5) 


The values at 500 °C, however, do not fall on the straight lines. This may 
be due to a low degree of accuracy in determining the Ml site occupancy par- 
ti cularly in Mg-rich samples at 500 °C. We may extrapolate and calculate the 
constants at 1000°C. An isotherm at 1000°C plotted with these constants is 
consistent for the most part with Virgo and Hafner’s (1969) data. 

Thermodynamic functions of mixing 

The 'partial 1 functions of mixing in Ml and M2 may be calculated by using 


the following relations: 

£M1 _ ^Ml xMl jj/Ml 
KM F ft wig 

(1.22) 


- = *“> a r M1 /a t 

(1.23) 
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similarly for the M2 site. 


The excess free energy of mixing for the crystal as a whole is given by 



The values of the activity coefficients / may be obtained by 


(7. 





(1 - x « 2)2 



x 


Opx 

Fe 


(7. 


similarly for f°p* . 


It is useful to represent the excess free energy of mixing as a polynomial 
in * as suggested by Guggenheim (1937): 



*sr> + 4 (*?r - *ar>* + > 


which may also be written as 


*Sr = A o + A i <*sr - ^ c*?r - *S5*) 2 + } <7. 


Equation (7.8) is of the form 


4 = a Q +. a l + a ? + 


and may be solved by a least squares program. To determine the constants A , 

A t and A 2 , we first obtain G°p x by using (7.6) and (7.7). The constants listed 
in Table 7.1 are used in (7.7). These values of excess free energy are then used 
to calculate by using a least squares program on equation (7.8) the three constants 
A q , A 1 and A 2 shown in Table 7.2. Figure 7.4 shows three curves at 600, 700 
and 800°C fitted to the data on reduced excess free energy {(7,4 * ) versus 

r e Mg 

(*F e " \ g )• This method of determining A Q , A 1 and A 2 is more reliable than 
the Redlich and Kister method followed by Saxena and Ghose (1971). 


The three constants as a function of l/T (Fig. 7.5) are given by 


A 0 = 10802 - 33862 


/i° 3 \ 

\ T 1 


+ 35135 ~ 11202 


(f)‘ 


(7.9) 


A 1 = 1789 - 3612 (I2?j + 2008 fly' 


(7.10) 


4 = - 13863 + 41051 - 40299 + 13426 (~) ( 7 - 1:L ) 


The heat of mixing and the entropy of mixing can now be readily calculated 
using the following equations: 


- S 0pX = *°P X r°P x 
°EM - *Fe X Wg 


B 4 ^ 4 
b r + b r 


OK X - *S') + 




*sp 2 + 


(7.12) 



(7.13) 
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The free energy of mixing is given by 


G °r = H %r - T s E°r + * T c*?r ^ + *sr ^ *sp (7 - M) 

It may be useful to write an expression for the free energy of mixing of the 
crystal as a whole in terms of partial thermodynamic functions of mixing 
referring to the sites. We have for (Mg, Fe) Si 2 O e accorciing to (6.23) 



-Ml - 

*Fe 


2 


x 


M2 

Fe 




(7.15) 


Note that 11°?* 4 (ffjg + ffg) and similarly 4 (Sg +S«2 ). The ex- 

pression (7.15) divides the total free energy of mixing into three terms. The 
term in the first bracket is the potential due to the difference in site occupancy 
energy. It includes both the contributions of the standard enthalpy and entropy 
of the exchange. The terms in the second and third brackets include the heat of 
mixing and the entropy of mixing respectively at the individual sites. 


This division of the free energy of mixing is artificial. It may be useful in 
understanding the change in free energy as a function of temperature and order- 
disorder. Figuers 7.6 and 7.7 show a comparison of the various functions of 
mixing at 600 and 800°C. The values of H m and TS m for the crystal as a whole 
are obtained by using equations (7.12) and (7,13). It may be noted that the first 
energy term 


x m _ * m 2 

Fe Fe 


A G° 

2 7a > 
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along with the entropy terms increase the negative C° px and make the solution 
stable. The negative value of the first energy term decreases somewhat with 
increasing temperature indicating that the disordering works against stability 
of the solution-. This effect is, however, limited and is overridden by the positive 
excess heat of mixing at the sites. 

At 800° C the contributions of the 'partial* energies of mixing at the sites 
to the energies of mixing in the crystal as a whole are somewhat asymmetric. 
The mixing energies lP px and 7 , S 0px have more or less symmetric values. 

EM EM 

With decreasing temperature, the energies of mixing at the sites increase, those 
for the Ml site increasing relatively more than those for the M2 site. The //° px 
and , therefore, become more and more asymmetric. 

Sites as Solutions With Quasi-chemical Approximation 
Equilibrium constant, K /M1 and IF M2 

As shown in section (5), the relation among K 7a , and the site 

occupancies may be expressed by the following equation: 



(7.16) 


where /3 and / 3 ' are for Ml and M2 respectively and given by: 

/3 = {1 + 4 0M1 0M1 ( e 2W m /zRT _ i^yi/2 
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P = {1 + 4 (e*" 2 /-® 1 - 1)}1/J 







4>'s are, as before (1.35), 

^Fe ^Fe 

*Fe ^Ffi + *Mg ?Mg 

z and z f are the number of neighboring cations in Ml and M2 respectively and 
q *s are contact factors. 

z and g for Ml and M2 sites 

As the silicate framework ^ jes not significantly change its character as a 
function of temperature or Fe 2+ /Mg 2+ ratio, the sites Ml and M2 have definite 
configurations and the polyhedral geometry changes only somewhat with changing 
Fe 2+ /Mg ratio. We may, therefore, assume that the number of sites which are 
neighbors of any one Ml site are two or four other Ml sites and q mi is equal 
to §mi and both are unity. A reference to Figure 7.1 shows that two of the 
four Ml sites are somewhat nearer to a central Ml. It may be noted that the 
two inner strips with Ml sites and two outer strips with M2 sites lie more or 
less in a plane. Since the two M2 strips are separated by the intervening Ml 
sites, it is only realistic to consider that number of neighboring M2 sites to any 
one M2 site is only two. We may also assume that <?“ 2 = ?“ 2 •- 1 and substitute 
in (7.16) z = z' = 2 and all q ’s as unity. 

Results of calculations 

The intracrystalline distribution data presented in Figure 7.2 was used to 
determine the quasi-chemical parameters W m and and the equilibrium 
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constant K 7a by using least squares analysis on equation (7.16). The distribution 
data at 500°C was not used and the number of data point at 700°C was not found 
sufficient for a satisfactory convergence. The values of H' M2 and K Ja at 600 
and 800°C are listed in Table 7.3. Chi square values for the simple mixture 
model and for the quasi-chemical approximation are not significantly different. 

Activity-composition relation at 600 °C 

We may calculate the ’partial’ activity coefficients at the Ml site by using 
the relations 


f Ml _ 
Fe 


<1 + 


< - !) \ 

(fi + i)i 


/2 


and similarly for f^ and finally determine the activity-composition relation in 
the crystal as a whole by the equation: 

4T = 0 $ f ") 1/2 ($1 f Fe ) 1/2 

The results of such calculations are presented in Figure 7.8. The figure 
shows that there is very little difference between the activities calculated by 
assuming sites as 'simple mixtures’ or as solutions with quasi-chemical 
approximation. 

Unmixing 

Although no data on unmixed orthopyroxenes are available, a possible un- 
mixing can be predicted theoretically. Due to kinetic reasons , data on actual 
unmixing may not be available. However orthopyroxenes of suitable compositions 


heated at temperatures close to the critical temperature of unmixing may show 
certain structural characters such as domains by the use of electron microscopy. 


The conditions of critical mixing are 

3 2 G E /3 x 2 = - R T/x (1 - x) 


(2.3) 


3 2 G E /B x 3 = - R T (2 x - l)/x 2 (1 - x ) 2 < 2 * 4 ) 

Successive differentiation of equation (7.10) leads to equations which are 
transcendental in character and cannot be solved without an iteration procedure. 

For the present purpose it is enough to find an approximate temperature where 
the activity-composition diagram shows the same activity over a small composition 
interval (Figure 7). From Figure 7.9 we find that the critical temperature is 
in the vicinity of 450 °C and the critical composition is about xOp* - 0.83. 

At 450°C, the values of H /M1 /RT and ff /M2 /RT are 2.23 and 1.49 respectively. 

For 0.83, we have = 0.75 and * ” 2 = 0.92. 

If Ml site was an independent phase, it would have unmixed into two coexist- 
ing phases with the following compositions and proportions approximately: 

Ml (Fe) 0.8 o'] Ml (Fe) 0.20*] 

>92% >8% 

Ml (Mg) 0.20J Ml (Mg) 0. 8()J 

For the above values of *mi x m remains nearly the same and we have 

Fe Fe 

Fe Ml -0.80 M2 0.9281 

>9 2% 

Mg Ml 0.20 M2 0.072 
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Fe Ml 0.20 M2 0.92 


(8% 

Mg Ml 0.80 M2 0.08 j 

The first of the above two represents the right direction towards unmixing. 
However, we cannot explain with above model, how a complete unmixing results 
at this critical temperature and composition. If the unmixing takes place via 
domain formation, it is obvious that the best experiments would be to examine 
the crystals with * 0px between 0.80 to 0.85 and which are heated between 450 

Fe 2 

to 500°C. 
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Table 7.1 


Values of Equilibrium Constant K 7a , AC° a and the Constants ff /M1 and 1P M2 


T°C 

*7 a 

w m 

Cal/mole 
W M2 

*G 7 ° a 

500 

.279 

2893 

1919 

1959 

600 

.277 

2390 

1577 

2233 

700 

.273 

1916 

1215 

2509 

800 

.289 

1641 

1057 

2646 

900 

.298 

1338 

834 

2818 

1000 

.311 

1102 

670 

2949 




Table 7.2 


The Three Constants A Q , A t and A 2 Obtained by a Least Squares Fit 
to Equation (7.8). The Values in Brackets are Values 
Calculated by Using (7.9-7.11) 


T°C 

'<0 

4 

a 2 

500 

1542 (1544) 

494 (478) 

871 (867) 

600 

1277 (1278) 

233 (287) 

454 (461) 

700 

976 (951) 

251 (199) 

321 (335) 

800 

652 (693) 

165 (168) 

298 (260) 

900 

548 (529) 

159 (170) 

143 (163) 
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Table 7.3 

Quasi-chemical Parameters ff' M1 and W m and K^ a as Calculated by (7.16) 


T°C 

K 7. 

IP M1 
R T 

ff /M2 

R T 

QC 

* 

Chi 2 

Simple 

mix. 

600 

0.293 

1.71 

1.26 

.027 

.017 

800 

0.273 

0.733 

0.469 

.007 

.004 


* Chi ^ " 2 ^ ( co ^ cu ^ g<,e< ^) ~ & (by least squares)} 2 

K (by least squares) 
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Figure 7-2. Distribution isotherms for the Mg y ' -Fe z partitioning between M 
in heated natural orthopyroxenes (Saxenc and Ghose, 1971). Th« 
least squares fit to the experimental data using ‘simple mixture' 
sites. 







Figure 7 •-'* ’’educed excess free energy of mixing plotted against (re p o — } in 

^rthopyroxene. The curves are leastsquares fit using the equation (7.8)- 
This fit yields the constants Aq, A , and 






Ca!/mo1e 



Figure 7*6. Thermodynamic functions of mixing in orthopyroxene at 600 °C, 
both at the sites and in the crystal as a whole plotted against 
#pP*. The functions of mixing at the sites correspond to the 
atomic fractions at Ml and M2 which in turn correspond to an 
equilibrium distribution at 600 °C in the crystal with x*?px as 
measured on the abscissa. A G ' is 
{(x Ml _ x ^ 2 )/2> A 










Section 8 


CRYSTALLINE SOLUTIONS AND GEOTHERMOMETRY 
The concept of metamorphic facies evolved through the attempts of the 
petrologlsts to distinguish the mineral assemblages formed at different P and 
T in the field. In several cases, the experimentalists have simulated the physical- 
chemical conditions of the formation of rocks in the laboratory. However, relev ent 
thermodynamic data on the rock forming minerals is meagre and causes difficulties 
in making a meaningful interpretation and in checking the validity of most experi- 
ments. From our existing knowledge, it is only possible to obtain certain quali- 
tative to semi-quantitative estimates of P and T of the formation of a mineral 
assemblage. Such methods are based on the knowledge of the chemical reactions 
which occur as a result of changing P and T within and in between the crystalline 
solutions. 

This section is not a review of the methods of geotherrnometry. Ingerson 
(1965) has made an excellent review of such methods. We shall be mainly con- 
sidered here with coexisting crystalline solutions and their thermodynamic 
properties of mixing as a function of P and? 1 . 

Inter-crystalline Equilibria 

Distribution of a component between coexisting simple mixtures 

Let us assume that we have data on the distribution of a component A between 
two binary phases (A, B) M and (A, B) N. The distribution data covers a large 
compositional range in both the phases. Let (A, B) M be a simple mixture whose 
activity-composition relations are known at several different temperatures but 
the activit 3 r -composition relation in (A, B) N is unknown. If the form of the 
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distribution isotherm, whose temperature is unknown, does not indicate (A, B) N 
to be too non-ideal, we may assume that the distribution follows the following 
equation: 


In K = In K„ + |1(1 - 2 *«.) - (1 - 2 *g.) < 8a ) 

where K is the equilibrium constant, the distribution coefficient and a and/3 
represent (A, B) M and (A, B) N respectively. 

To find the temperature of the distribution isotherm (considering that the 
effect of P is negligible on the equilibrium constant K ), we may assume the 
temperature to be T and solve equation (8.1) by a least squares method. We 
repeat this work by assuming various temperatures T , T etc. In practice we 
could rely on geologic evidence to fix the limits T l and T 2 within which the 
actual temperature of the ion- exchange equilibrium might lie. These calculations 
yield a set of W a and values at temperatures 7’ 1 and T . Such values of W a 
at T 1 and T 2 differ from the experimental values of w a at these temperatures 
because of the actual temperature of the isotherm being different from the 
assumed T l or T 2 . We may, therefore, plot the calculated values of W a and the 
experimental values of W a against T as in Figure 8.1. The intersection of the 
two lines yields the temperature of the distribution isotherm. 

This method was used by Saxena (1971) in estimating temperature of crys- 
tallization for coexisting pyroxenes in the metamorphic rocks. Unfortunately 
most minerals are not binary solutions and many of those which are quasi-binary 
may not be strictly symmetric simple mixtures. 


Distribution of NaAlSi 3 O g between coexisting plagioclase a nd alkali feldspar 
Barth (1956, 1962) proposed that the ratio mole fraction of albite (Ab) in 
alkali feldspar: mole fraction of albite (Ab) in plagioclase varies linearly with 
inverse of absolute temperature. This is true if mixing of the end members 
Ab and Or (orthoclase, KAlSi 3 O g ) and Ab and An (anorthrite, Ca 2 Al 2 Si 2 O g ) is 
ideal or nearly ideal in the compositional range. 


The properties of sanidine crystalline solution has been studies by Thompson 
and Waldbaum (1968, (1969) and by Thompson (1969a). Even in the high temperatu: 
range 700 to 1000°C, the solution deviates significantly from ideal mixture. 
Similarly it is expected that there will be significant deviation from ideality 
in the plagioclase feldspars. 


Plagioclase and sanidine may be Considered as two ternary solutions at a 
certain P and T, we have the following ion-exchange reactions 


agioclase Sanidine 
NaAlSi 3 O g ^NaAlSi 3 O g 

(8.2) 

K AlSi 3 O g AlSi 3 0 8 

(8.3) 

CaAl 2 Si 2 O g ^ CaAI 2 Si 2 0 g 

(8.4) 


The corresponding equilibrium constants are 


K n = 


„s f s 

-S/-B =— — 


,P1 f Pi 
I -A 


Ab ' Ab 


fU = 


*8,. 


a "_/ 1 


Pi 

Or 


X n fn 
Or 1 Or 

* P1 f pl 

Or 1 Or 
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k a 


s 

An 


f s 

‘ An 


— ■■ • 


pi 

* P1 

/■pi 

An 

*An 

' An 


where x = Ab/Ab + An + Or. The activity coefficients f are functions of P , T and 
the ratio of two of the other mole fractions. 


In actual practice, we may simplify the situation somewhat by assuming 
that both sanidine and plagioclase are binary solutions. Therefore we write K x 
as 


K 


2 




/lb 


where *® = Ab/Ab + Or in sanidine and - Ab/Ab + An. ff, is now a function 
of P , T and *| b , and similarly / p * . Even in the simplified case, we have to know 
the activity-composition relations in both sanidine and plagioclase solutions. 

The activity-composition data on sanidine crystalline solution is available 
(Thompson and Waldbaum 1968, 1969b) but we have yet to learn about the pla- 
gioclase crystalline solution. Only sanidine and high temperature plagioclase 
have been considered above. The temperature estimate at lower range of 
temperatures is further complicated by the structural changes in the feldspars 
as a function of T and composition. 


Distribution of Fe and Mg in Coexisting Minerals 
The use of the distribution coefficient involving two multicomponent minerals 
in problems of petrogenesis is severely limited. The distribution of Fe 2+ and 
Mg 2+ between coexisting garnet and biotite is one example. Such a distribution 
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is a function of P , T and the concentrations of Mn and Ca in garnet and of Al 3+ , 
Fe 3+ and Ti in biotite (Kretz, 1959, Albee, 1985, Sen and Chakraborty, 1968, 
Saxena 1968 c). 

Fortunately garnet-biotite pair is one of the most investigated pairs and 
several chemical analyses are available in literature. It is, therefore, possible 
to improve the usability of the by analyzing its compositional dependence by 
a multivariate statistical analysis. Saxens (1969b) used the ’principal component 
analysis' and obtained the ’transformed distribution coefficient’ as 

Transf. = 0.5013 - 0.4420 xf° r + 0.1506 - 0.3474 x*» r 

+ 0.0865 x** 1 ~ 0.0333 * bi TV - 0.3165 * bi VT 

Ca A1 IV Ai VI 

+ 0.5488 * b [ 

where * pe is Fe 2+ /(Fe 2+ + Mg), K D is x b * x*>® r /x|® r xjj and the other *’s rep- 
resent ions on 12 O basis in garnet and 22 O basis in biotite. The transformed 
Kq could be shown to vary systematically with temperature inferred on geologic 
evidence. Figure 8.2 shows the use of such a transformed in estimating the 
temperature of crystallization of rocks containing garnet and biotite. Using 
this method of estimating temperature, Dennen, Blackburn and Quesada (1970) 
obtained an average temperature of 650 °C for certain Grenville gneisses which 
agrees well with their independent estimate based on the concentration of Al 3+ 
in quartz. Some of the garnet compositions particularly from low grade met- 
amorphic rocks may represent zoning. Therefore, there is further possibility 
of improving the quality of the transformed K D as a temperature indicator. It 
is assumed that P variation between 4 to 8 Kb does not change and other 
substitutional relations .significantly. 
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Concentration of Ca 2+ in clinopyroxene 

The concentration of Ca 2+ in clinopyroxene coexisting in equilibrium with 
orthepyroxene may be used to estimate the P and T of crystallizstion of the 
rock. The solvus relations in the systems enstatite diopside (Davis and Boyd, 

1966, Warner, 1971) and ferrosilite-hedenbergite (Lindsley and Munoz, 1969) may 
be used for this purpose. Change in P has only slight influence on this relation- 
ship (Davis and Boyd, 1966a, Warner, 1971). The main uncertainty in such esti- 
mates arise due to the fact that climopyroxene in natural assemblages are not 
binary solutions. Even if the concentrations of Al 3+ , Fe 3+ and Ti 4+ are neglected, 
we still have to reckon with a ternary solution with Fe 2+ , Ca 2+ and Mg 2+ , 

Intra-crystalline Equilibria 

Order -disorder 

Chemical reactions within the crystal involve the order-disorder phenomena 
or the distribution of cations in the nonequivalent structural sites. The distri- 
bution isotherms for Fe 2+ -Mg distribution between Ml and M2 sites in ortho- 
pyroxene shown in Figure 7.1 may be used to estimate the temperature of 
equilibrium for the intracryst'alline ion-exchange reaction. Figure 8.3 shows 
a geothermometer for this purpose. A knowledge of the distribution coefficient 
for the ion-exchange equilibrium 

Fe 2+ (M2) + Mg 2+ (Ml)^ Fe 2+ (Ml) + Mg 2+ (M2) 
and x 2 * in the orthopyroxene is enough for estimating the equilibrium temperature. 

The intracrystalline ion-exchange is a rapid process. The ion-exchange, 
therefore, continues even after the crystallization has been completed and until 
the rock has cooled below a certain critical temperature. If the rock is quenched 
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as certain meteor itic bodies and some volcanic lavas, the intracrystalline 
equilibrium temperature would correspond to the temperature of the formation 
of the rock. In all other cases this temperature will be much lower than the 
actual temperature of crystallization. More information on the cooling rates of 
the rocks and the kinetics of ion-exchange process is necessary for estimating 
the temperature of crystallization of the mineral assemblage. 

The peristerite solvus 

Many crystalline solutions show unmixing below a certain critical temperature. 
This property may sometimes be used as an indicator of P and T. Plagioclases 
(An 2 to An 18 ) unmix into two phases with comppsitions in the range An 0 _ 1 and 
An 2 s - 28 a certain critical temperature. This temperature has not been 

determined but it probably lies below 600 °C which is the temperature of Al-Si 
ordering in the plagioclase (see Ribbe I960, 1962, Brown, 1960, 1962). An ex- 
perimental determination of peristerite solvus could be very useful for estimating 
temperature of crystallization. Crawford (1966) has made a useful study of the 
composition of plagioclases in rocks of varying metamorphic grade. Siich pet- 
rologic studies should be useful to those who are attempting to fix the peristerite 
solvus experimentally. 

In the system FeS-znS, the amount of FeS in sphalerite is a function of T of 
formation (Kullerud, 1953) when FeS coexists as a phase with sphalerite. 

Sub-solidus compositional relationships between coexisting phases such as 
calsite-dolomite (Goldsmith and Newton, 1969) and forsterite-monticelite (Warner, 
1971) may also be used for temperature estimates. It should be possible to 
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consider unmixing in some other crystalline solutions such as magnetic (see 
Rumble, 1970) as probable indicators of P and T in natural assemblages. 
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